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ABSTRACT 


7&The  purpose  of  this  study  was  twofold:  to  evaluate  the  ability  of  a 
circularly  symmetric  two-cylinder  model  of  a precipitating  cumulus  cell 
to  simulate  radar-observed  cells,  and  to  ascertain  to  what  accuracy  the 
model  input  variables  need  to  be  known  to  obtain  realistic  results.  ■£< 

Data ''from  GATE  were  employed  because  of  their  availability  and  cover- 
age andfbecause  the  tropical  atmosphere  during  the  experiment  contained 
relatively  little  vertical  wind  shear.  Cells  organized  into  two  types 
of  mesoscale  features,  a line  and  a cluster,  were  simulated.  Because  of 
the  variability  observed  in  the  cell  characteristics,  the  mesoscale 
features  were  idealized  into  geometrically  simple  configurations  of 
'representative*  cells.  The  representative  cells  were  geometrically  com- 
patible with  the  model  and  had  dimensions  based  on  the  average  sizes  of 
observed  cells.  The  size  of  the  model  inner  cylinder,  assumed  to  be  the 
region  of  active  convection,  was  estimated  from  the  areas  of  observed 
reflectivity  cores.  The  size  of  the  outer  cylinder,  analogous  to  the 
environment  of  the  cells,  was  estimated  from  distances  between  cells  in 
the  mesoscale  features.  This  idealization  process  produced  what  appeared 
to  be  realistic  dimensions. 

Vertical  profiles  of  temperature  and  relative  humidity  were  also 
needed  for  initialization.  A large  sample  of  observed  profiles^was  ideal- 
ized by  a method  which  retained  the  main  characteristics  of  the  profiles 
while  smoothing  the  minor  variations.  The  idealized  relative  humidity 
profiles  fell  into  three  distinct  types,  the  distribution  of  which  varied 
with  the  synoptic  situation.  Variations  in  the  idealized  temperature 
profiles  were  small  and  not  systematic. 

The  effectiveness  of  the  model  in  simulating  observed  cells  was  in- 
vestigated by  comparing  observed  with  model  cell  heights,  maximum  reflec- 
tivities, and  heights  of  the  reflectivity  maxima.  When  initialized  appro- 
priately, the  model  produced  cells  which  appeared  to  be  remarkably  similar 
to  observed  cells.  Although  cell  dimensions  were  not  unimportant,  the  rela- 
tive humidity  profile  exerted  the  greatest  influence  on  model  cell  growth. 


P A ^ P 
t-r  n }P 


- 3 - 


Open-atmosphere  humidity  profiles  were  found  to  be  too  dry  to  support 
vigorous  convection.  Realistic  cells  were  generated  when  the  humidity 
was  increased  under  the  assumption  that  the  cells  were  growing  in  a 
mesoscale  region  which  was  more  moist  than  the  open  atmosphere.  Calcu- 
lations showed  that  convective  activity  could  contribute  to  the  humidi- 
fication of  the  open  atmosphere,  but  this  contribution  was  not  sufficient 
to  alter  the  relative  humidity  profile  to  the  point  where  realistic  cells 
could  be  simulated.  Further  humidification  via  mesoscale  ascent  was 
necessary. 

The  need  for  taking  into  account  the  effects  of  mesoscale  processes 
is  discussed,  and  methods  for  doing  so  are  suggested.  The  need  for  the 
addition  of  a realistic  boundary  layer  is  discussed.  The  further  addi- 
tion of  a less  artificial  initiating  impulse  coupled  with  the  structure 
and  properties  of  the  boundary  layer  is  suggested. 
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1.1  Background 

It  is  generally  recognized  that  moist  convection  plays  an  impor- 
tant role  in  vertical  heat  and  moisture  transport.  Through  this  role,  it 
is  believed  to  have  far-reaching  effects  on  synoptic  and  global  circula- 
tions (see,  e.g.  Charney  and  Eliassen,  1964;  Gray,  1973).  These  effects 
cannot  be  completely  understood,  however,  until  the  nature  and  magnitude 
of  cumulus  transports  can  be  determined  with  a reasonable  degree  of  accur- 
acy. Because  of  the  small  scale  of  cumulus  cells  and  their  complexity, 
it  is  extremely  difficult  to  measure  directly  their  transport  of  heat, 
moisture,  and  other  qualities.  Indirect  methods  of  determining  cumulus 
transports  using  radar  data  (Austin  and  Houze,  1973)  or  diagnostic  studies 
of  large-scale  variables  (Arakawa  and  Schubert,  1974)  have  been  devised, 
but  they  provide  only  crude  estimates. 

One  tool  which  can  aid  the  meteorologist  in  investigating  cumulus 
transports  is  the  cumulus  cell  model.  Cell  models  have  been  designed  to 
serve  a number  of  purposes  (investigation  of  microphysical  or  dynamical 
processes,  aids  in  weather  modification  experiments,  etc.).  Some  of  these 
models  might  be  suitable  for  use  in  cumulus  transport  studies  provided 
that  they  are  capable  of  simulating  observed  cells  with  a sufficient 
degree  of  detail.  Comparisons  of  cells  generated  by  the  more  theoretically 
oriented  models  with  observed  cells  have  usually  been  made  in  the  most 
general  way,  involving  only  the  concern  that  the  model  cell  resemble  a 
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(l.e.  any)  cumulus  cell  (see,  e.g.  Orville  and  Sloan,  1970,  who  considered 
the  configuration  of  orographically  induced  cells).  In  other  studies 
(e.g.  Weinstein,  1970),  in  which  particular  observed  cells  were  modeled, 
qualitative  comparisons  have  been  made  between  the  general  appearance  of 
major  features  in  the  modeled  and  observed  cells.  The  quantitative  com- 
parisons that  have  been  made  (e.g.  Sax,  1969,  who  compared  computed  with 
observed  cell  heights)  were  not  made  in  sufficient  detail  to  elicit  confi- 
dence in  the  models'  ability  to  simulate  important  cumulus  processes. 

The  main  difficulty  in  making  comparisons  has  stemmed  from  the 
lack  of  important  processes  or  properties  which  could  be  both  observed  in 
a real  cell  and  calculated  in  a model  with  sufficient  accuracy  and  reso- 
lution to  allow  a reasonably  accurate  comparison.  Yau  (1977)  has  contri- 
buted to  the  solution  of  this  problem  by  devising  a microphysical  para- 
meterization which  allows  for  the  natural  development  of  drop  size  spectra 
and  thus  the  calculation  of  realistic  radar  reflectivities.  These  reflec- 
tivities can  be  compared  with  reflectivities  observed  within  real  cells 
by  radar.  Yau  incorporated  his  microphysical  parameterization  into  a 
fairly  sophisticated  two-cylinder  model  in  which  an  inner  cylinder  repre- 
sents a cumulus  cell  and  a surrounding  outer  cylinder  represents  the 
environment  of  the  cell.  This  type  of  formulation,  while  being  reasonably 
well  suited  for  investigating  cumulus  transports,  requires  less  computation 
than  the  more  complex  two-dimensional  models  of,  for  example,  Soong  and 
Ogura  (1976).  This  computational  simplicity  permits  relatively  extensive 
experimentation  to  be  carried  out. 
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1.2  Statement  of  Problem 

In  spite  of  the  advantages  of  Yau's  model,  it  is,  compared  to  the 
real  world,  highly  simplified.  In  order  to  assess  its  usefulness  for  com- 
puting vertical  transports,  it  is  desirable  to  determine  how  well  it  simu- 
lates real  cells.  Toward  this  objective,  this  study  will  comoare  cells 
generated  by  the  model  with  cells  observed  by  radar.  A secondary  objective 
is  to  determine  specifically  what  characteristics  of  cells  and  the  atmo- 
sphere  in  which  they  are  occurring  need  to  be  known,  and  to  what  accuracy 
they  need  to  be  known,  in  order  to  achieve  realistic  simulations  from  the 
model . 

The  model  will  be  initialized  with  horizontal  cell  dimensions 
estimated  from  radar  data  and  with  environmental  temperature  and  relative 
humidity  profiles  determined  from  rawinsonde  soundings.  Comparisons  will 
be  made  between  radar  reflectivities  generated  by  the  model  and  those  ob- 
served by  radar.  The  data  will  be  drawn  from  the  GARP*  Atlantic  Tropical 
Experiment  (GATE).  This  data  set  was  chosen  for  its  availability  and  ex- 
cellent coverage  and  because  the  tropical  atmosphere  has  relatively  little 
vertical  wind  shear,  a property  which  cannot  be  taken  into  account  by  the 
axisymmetric  cell  model. 

In  Chapter  2,  the  model  microphysics,  dynamics,  and  general  formu- 
lation are  briefly  described.  The  input  variables  are  discussed  and  the 
results  of  Yau's  sensitivity  tests  on  those  variables  are  presented. 

Values  of  input  variables  which  will  remain  constant  throughout  this  study 
and  the  rationale  behind  their  choice  will  be  given. 
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The  data  are  discussed  in  Chapter  3.  A brief  description  of  GATE 
is  given  after  which  the  sources,  handling,  and  analysis  of  the  data  are 
described.  The  need  for  idealizing  the  data  to  provide  "representative" 
cell  dimensions  and  soundings  is  discussed  and  the  idealization  process 
is  described.  Methods  for  comparison  of  model -generated  with  radar- 
observed  cells  are  presented. 

Two  types  of  situations  are  considered.  One  is  a well -organized 
and  well -developed  line  of  convective  cells,  a configuration  which  is 
typical  of  tropical  disturbances.  Experiments  and  comparisons  for  this 
situation  are  presented  in  Chapter  4,  The  second  type  of  situation,  a 
cluster  of  cells  which  occurred  during  a relatively  undisturbed  period, 
is  considered  in  Chapter  5. 

Conclusions  drawn  from  the  experiments  and  recommendations  for 
future  use  of  the  model  are  presented  in  Chapter  6, 


CHAPTER  TWO 


DESCRIPTION  OF  MODEL 

The  model  being  evaluated  in  this  thesis  was  designed  by  M.K.  Yau 
(1977)  to  allow  investigation  of  the  interaction  of  dynamical  and  micro- 
physical processes  in  precipitating  cumulus  clouds.  In  this  chapter,  a 
brief  description  of  the  important  characteristics  and  properties  of  the 
model  is  given. 

2.1  Microphysical  Processes 

The  microphysical  processes  were  designed  to  be  a compromise 
between  the  explicit  calculation  of  diffusional  growth  of  drops  and  of 
stochastic  collection  processes  and  the  simplified  parameterization 
scheme  proposed  by  Kessler  (1969),  in  which  interaction  between  cloud  and 
rain  is  formulated  on  the  concept  of  autoconversion  and  in  which  precipi- 
tation particles  always  assume  the  Marshall -Palmer  distribution  and  fall 
within  a single  effective  velocity  at  a given  level.  As  in  the  Kessler 
parameterization,  Yau  bypassed  growth  of  cloud  droplets  and  diffusional 
growth  of  precipitation  particles  by  employing  an  autoconversion  scheme 
and  allowing  for  accretion  of  cloud  water  onto  precipitation  particles. 
However,  he  grouped  precipitation  particles  into  ten  size  categories, 
each  with  a representative  fall  velocity,  and  required  that  autoconverted 
water  be  added  only  to  the  smallest  size  category.  The  drop-size  distri- 
bution was  then  built  through  accretion  of  cloud  and  interaction  between 
precipitation  particles.  The  efficiency  for  the  collection  of  cloud  drop 
lets  by  precipitation  particles  was  assumed  to  by  0.8.  The  collection 
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efficiency  for  raindrops  was  based  on  a collision  efficiency  of  one  and  the 
coalescence  efficiency  given  by  Brazier-Smith  j?t  a2.  (1973).  Breakup  pro- 
cesses included  were  impaction  (in  which  each  collision  produces  a main 
drop  and  three  satellite  drops)  and  instantaneous  breakup  (in  which  drops 
larger  than  those  in  the  largest  size  category  spontaneously  break  into  a 
number  of  drops  distributed  over  a pre-determined  spectrum).  Finally, 
cloud  existing  in  subsaturated  air  was  assumed  to  evaporate  instantaneous- 
ly to  the  point  of  saturation.  If  the  cloud  water  was  insufficient  to 
saturate  the  air,  rain  water  was  evaporated  using  the  microphysical  dif- 
fusion equation. 

Yau  tested  his  microphysical  parameterization  against  the  stochas- 
tic scheme,  comparing  the  initiation  of  warm  rain  in  a one-dimensional 
updraft.  Comparisions  between  liquid  water  content  and  radar  reflectivity 
at  two  heights  in  the  updraft  and  of  surface  rainfall  intensity  showed  no 
significant  differences  between  the  two  schemes.  Comparisons  of  the  evo- 
lution of  raindrop  size  spectra  showed  a few  minor  differences  of  question- 
able significance. 

2.2  Model  Formulation  and  Dynamics 

The  model  was  designed  to  employ  the  anelastic  pressure  equation 
(Ogura  and  Phillips,  1962)  in  such  a way  as  to  allow  explicit  calculation 
of  the  perturbation  pressure.  Geometrically,  the  model  is  composed  of  two 
concentric  cylinders  whose  radii  are  constant  with  height  and  time.  The 
inner  column  is  the  cloud  region  which  comprises  the  area  of  active  con- 
vection. Compensatory  currents  occur  in  the  outer  column  to  satisfy  the 
requirement  of  continuity.  Horizontally  averaged  equations  were  developed 
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in  the  manner  described  by  Asai  and  Kasahara  (1967),  in  which  the  cloud 
is  assumed  to  be  axi symmetric  and  the  properties  at  a given  level  are 
averaged  across  each  cylinder.  Exchange  of  properties  between  cylinders 
is  accomplished  by  dynamic  entrainment  and  detrainment  (i.e.  one-way 
advection  by  horizontal  flow  between  cylinders)  and  turbulent  entrainment, 
a two-way  exchange  of  properties,  which  is  a function  of  the  difference 
of  those  properties  and  the  difference  of  the  vertical  motion  across  the 
interface.  Aside  from  this  type  of  turbulent  exchange,  vertical  and  hori- 
zontal eddy  fluxes  are  neglected.  Boundary  conditions  specify  no  flow 
through  the  top  or  bottom  of  the  cylinders  or  through  the  outer  vertical 
boundary  of  the  outer  cylinder.  The  outer  cylinder  thus  becomes  the 
"environment"  and  is  isolated  from  interactions  with  larger  scale  effects. 

The  model  uses  specifications  of  initial  temperature  and  relative 

humidity  versus  height  to  define  a pressure-versus-height  profile  which 

results  in  a hydrostatic  atmosphere.  Convection  is  then  initiated  by 

imposing  a vertical  velocity  and  moisture  impulse  on  the  bottom  of  the 

2 

inner  cylinder.  The  velocity  impulse  takes  the  form  of  a sin  function  of 
height  while  the  moisture  impulse  saturates  part  or  all  of  the  region  of 
velocity  impulse. 

2.3  Model  Input  Variables 

The  model  allows  the  specification  of  a number  of  variables  of 
meteorological  and  numerical  significance  in  order  to  simulate  cells  under 
different  conditions  and  handle  the  calculations  in  an  appropriate  way. 

The  meteorological  variables  specify  the  size  of  the  cell,  the  character- 
istics of  the  atmosphere  in  which  it  is  occurring,  the  nature  of  the 
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initiating  impulse,  and  a microphysical  property  related  to  the  conversion 
of  cloud  water  to  precipitation.  The  numerical  variables  deal  with  the 
numerical  space  and  time  resolution,  the  precision  of  calculations,  and 
the  time  limit  of  integration.  In  this  section,  the  variables  will  be 
described,  their  significance  will  be  discussed,  the  results  of  Yau's 
sensitivity  tests  on  them  will  be  summarized,  and,  when  possible,  the 
values  to  be  used  in  this  study  and  the  reason  for  their  choice  will  be 
given. 

2.3.1  Meteorological  Variables 

The  first  set  of  meteorological  variables,  specifying  the  size  of 
the  cell  to  be  modeled  and  the  size  of  its  environment,  includes  the 
radius  of  the  inner  cylinder  and  the  ratio  of  the  inner  to  outer  cylinder 
radii.  The  radius  of  the  inner  cylinder  defines  the  area  of  the  initia- 
ting updraft  and  hence  the  region  of  active  convection.  The  ratio  of  the 
inner  to  outer  cylinder  radii  (designated  o)  defines  the  radius  of  the 
outer  cylinder  given  the  radius  of  the  inner  cylinder.  According  to  the 
dynamics  of  the  model,  changes  in  the  inner  cylinder  radius  (without 
changing  o)  affect  the  strength  of  the  perturbation  pressure  force  and  the 
entrainment  process.  Both  the  dynamic  and  turbulent  entrainment  bear  an 
Inverse  relation  to  the  inner  radius.  The  perturbation  pressure  force,  on 
the  other  hand,  increases  as  the  radius  of  the  inner  cylinder  increases. 
Yau  conducted  experiments  using  inner  cylinder  radii  of  1.0,  1.5,  and  2.0 
km.  The  vigor  of  the  cell  increased  dramatically  as  the  radius  was  in- 
creased from  1.0  to  1.5  km,  but  more  slowly  as  the  radius  was  further 
Increased  to  2.0  km.  These  results  seem  to  indicate  that  entrainment  is 
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an  important  dissipative  process  for  the  smaller  cells  while  the  pressure 
force  takes  on  more  importance  in  the  larger  cells.  The  ratio  of  the 
radii  defines  the  ratio  of  the  areas  of  the  cylinders.  This  ratio  plays 
an  important  role  in  the  averaging  of  properties  over  the  cylinders  and 
in  continuity  considerations.  Proceeding  on  the  hypothesis  of  Asai  and 
Kasahara  (1967)  that  convection  will  maximize  vertical  heat  transport, 

Yau  calculated  vertical  heat  transports  for  reasonable  combinations  of  the 
inner  cylinder  radius  and  the  lapse  rate  while  varying  a.  The  results 
indicated  that  o should  be  between  0.3  and  0.4  (depending  on  the  situation) 
to  maximize  vertical  heat  transport.  In  all  his  experiments  except  one, 

Yau  used  a a of  0.32.  The  test  case  used  a o of  0.16  and  yielded  a much 
stronger  cell.  The  reasons  for  this  are  as  follows.  Continuity  requires 
that  the  updraft  that  occurs  in  the  inner  cylinder  during  the  development 
of  the  cell  be  compensated  by  a downdraft  in  the  outer  cylinder. 

Because  the  area  of  the  outer  cylinder  is  smaller  for  the  larger  a,  the 
compensatory  downdraft  must  be  stronger.  This  affects  the  growth  of  the 
cell  in  two  ways.  First,  the  mixing  of  downward  momentum  by  entrainment 
reduces  the  strength  of  the  updraft  and  hence  the  rate  of  condensation. 
Secondly,  adiabatic  warming  and  drying  of  the  air  in  the  outer  cylinder 
and  its  subsequent  mixing  inside  the  cloud  reduces  the  amount  of  vapor 
that  can  be  condensed.  The  result  is  a weaker  updraft  and  less  condensa- 
tion for  a larger  o.  In  this  study,  the  values  of  the  inner  cylinder 
radius  and  o will  be  based  on  radar  data. 

The  second  set  of  meteorological  variables  describes  the  atmosphere 
in  which  the  convection  is  taking  place.  The  height  of  the  model  is  speci- 
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fied  in  such  a way  that  the  top  of  the  model,  where  vertical  motion  is 
constrained  to  cease,  provides  as  little  interference  as  possible  with  the 
growth  of  the  cell.  This  depends,  of  course,  on  the  nature  of  the  cell 
being  modeled,  but  for  cells  which  are  not  expected  to  penetrate  the 
tropopause  it  is  reasonable  to  place  the  top  of  the  model  near  the  tropo- 
pause. Yau  used  15  km  for  the  model  height  in  his  tests  and  did  not  test 

the  model's  sensitivity  to  this  variable  (probably  because  none  of  his 
cells  reached  over  10  km  in  height).  This  study  will  also  use  15  km  for 
the  model  height  since  that  is  the  approximate  height  of  the  tropopause 

and  since  no  cells  were  observed  to  reach  that  height  in  the  situations 

being  tested. 

The  model  allows  the  assumption  of  an  idealized  atmosphere  or  the 
input  of  an  observed  atmosphere.  For  an  idealized  atmosphere,  the  vertical 
profiles  of  temperature  and  relative  humidity  are  constructed  from  given 
lapse  rates.  Describing  an  actual  atmosphere  requires  designating  a sur- 
face pressure  and  temperature  and  relative  humidities  at  500  m increments 
from  the  surface  to  15  km.  The  sensitivity  of  the  model  to  these  variables 
is  rather  obvious  — stronger  lapse  rates  and  higher  relative  humidities 
result  in  more  vigorous  cells.  In  this  study,  the  surface  pressure  will 
be  determined  from  nearby  surface  observations  and  the  vertical  profiles 
of  temperature  and  relative  humidity  will  be  determined  from  observed 
soundings. 

The  third  set  of  meteorological  variables  describes  the  initiating 
impulse.  These  variables  are  the  height  of  the  impulse,  the  strength  of 
the  velocity  impulse,  and  the  degree  of  saturation  of  the  impulse.  The 
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velocity  profile  in  the  impulse  is  given  by  w = wQ  sin^  (pp)  for  Z < Hq, 
where  HQ  is  the  height  of  the  impulse  and  wQ  is  its  maximum.  Yau  used 
an  Hq  of  2 km  in  all  his  experiments  and  got  reasonable  results.  Since 
this  impulse  scheme  is  rather  artificial,  no  attempt  will  be  made  to 
determine  a mors  appropriate  height.  Yau  used  a wQ  of  1 m sec”1  in  all 
his  experiments  except  one,  in  which  0.2  m sec”1  was  used.  The  resulting 
cell  in  the  test  case  was  nearly  identical  to  the  ’'control"  storm  except 
that  it  took  slightly  longer  (about  5 minutes)  to  develop.  Since  timing 
is  not  critical  in  this  study  (indeed,  the  time  interval  at  which  radar 
data  were  available  makes  such  a time  lag  insignificant),  a wQ  of  lmsec”1 
will  be  used.  The  impulse  may  be  saturated  over  its  entire  height  or  only 
above  the  cloud  base.  Yau  did  not  test  these  alternatives  in  his  thesis. 
He  did,  however,  test  a cell  saturated  above  the  cloud  base  against  a cell 
with  no  humidity  impulse.  The  development  of  the  latter  cell  was  delayed, 
but  otherwise  it  did  not  differ  significantly  from  the  cell  saturated 
above  the  cloud  base.  Since  timing  is  not  important,  this  study  will  use 
an  impulse  which  is  saturated  over  its  entire  depth. 

The  last  meteorological  variable  is  the  autoconversion  threshold. 
When  the  amount  of  cloud  water  exceeds  this  value,  autoconversion  takes 
place.  Yau  tested  autoconversion  thresholds  for  cloud  water  mixing 
ratios  between  0.0  and  0.001  in  a kinematic  updraft  and  compared  rainfall 
intensity  as  a function  of  time.  This  test  demonstrated  that  precipita- 
tion was  relatively  insensitive  to  changes  in  autoconversion  threshold. 

For  simplicity,  Yau  used  0.001  9 g"1  in  his  study.  The  same  value  will  be 
used  here. 


2.3.2  Numerical  Variables 

The  first  numerical  variable  prescribes  the  vertical  resolution 
of  the  model.  The  gridpoints  in  a 15  km  high  model  may  be  500  m or  250  m 
apart.  Yau  used  250  m in  all  his  experiments  except  one,  in  which  he  used 
500  m.  This  test  showed  that  the  larger  grid  spacing  slightly  decreased 
the  amplitude  of  most  cloud  variables.  Since  the  decrease  was  not  signi- 
ficant, a grid  spacing  of  500  m will  be  used  in  this  study  to  save  compu- 
tation time. 

The  second  numerical  variable  prescribes  the  time  step.  Yau  used 
10  seconds  in  all  his  experiments  except  two,  in  which  he  tested  5 sec  and 
20  sec.  The  results  demonstrated  that  the  model  is  not  sensitive  to 
changes  in  the  time  step  within  this  range.  A time  step  of  20  sec  will  be 
used  here  to  save  computation  time. 

The  third  numerical  variable  prescribes  the  number  of  Fourier 
modes  used  in  calculating  the  perturbation  pressure.  Yau  used  ten  modes 
in  all  experiments  except  one,  in  which  he  tested  twenty  modes  and  found 
no  significant  change.  Again,  to  save  time,  this  study  will  use  ten  modes. 

The  last  numerical  variable  prescribes  the  time  of  integration. 

This  should  be  long  enough  to  allow  the  cell  to  complete  most  of  its  life 
cycle,  but  not  so  long  as  to  waste  computation  time.  Experience  from  Yau's 
experiments  and  general  knowledge  of  thunderstorms  indicate  that  an  hour 
should  be  a sufficiently  long  time  to  integrate  the  model. 

2.4  Model  Output  Variables 

A myriad  of  quantities  related  to  dynamic,  thermodynamic,  and 
microphysical  processes  may  be  output  for  each  model  level  at  given  times 


i : 


during  the  course  of  the  integration.  The  output  variables  used  for  com- 


paring model  cells  to  observed  cells  will  be  discussed  in  Chapter  Three. 
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CHAPTER  THREE 

SOURCES  OF  DATA  AND  METHODS  OF  ANALYSIS 


3.1  The  Source  of  the  Data 


In  order  to  make  the  type  of  comparison  attempted  in  this  thesis, 
a reasonably  complete  data  set  — composed  of  surface,  upper  air,  and  radar 
observations  --  was  needed.  Routine  operational  data  would  not  have  pro- 
vided the  needed  detail;  a comprehensive  data  set  from  a special  research 
project  was  required.  Such  a data  set  was  available  from  the  GARP  Atlantic 
Tropical  Experiment  (GATE).  GATE  was  conducted  in  the  equatorial  region 
from  West  Africa  to  the  western  Atlantic  Ocean  during  the  summer  of  1974. 

To  meet  the  general  objective  of  investigating  the  mechanisms  by  whicn 
tropical  heating  drives  the  global  circulation  of  the  atmosphere,  70  nations 
participated,  providing  40  ships,  13  aircraft,  and  several  satellites  to 

serve  as  observational  platforms.  Meteorological  phenomena  were  investi- 

* 

gated  on  a number  of  scales  from  synoptic  to  cumulus.  The  scale  of  great- 
est interest  in  this  thesis  is  the  "cloud-cluster"  scale,  referred  to  as 

2 3 

the  B-scale,  which  encompassed  a range  of  10  to  10  km.  Data  for  inves- 
tigating phenomena  of  the  cloud-cluster  scale  were  gathered  primarily  by 
an  array  of  seven  ships,  referred  to  as  the  B-scale  array.  In  this  array, 
centered  at  8.5°N,  23.5°W,  the  seven  ships  formed  a hexagon  with  one  ship 
in  the  center  (see  Figure  3,1),  A variety  of  surface  and  oceanic  data 
was  recorded  at  each  of  the  ships  and  sondes  were  released  every  three 

1 

hours  (1-1/2  hours  in  especially  interesting  situations).  Four  of  the 


ships  carried  quantitative  weather  radars  with  digital  recording  systems. 
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Figure  3. 1.  The  GATE  ship  arrays  during  Phase  III.  The  inner  hexagon 
is  the  B- scale  array. 
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The  experiment  was  composed  of  three  21-day  observational  phases,  Data 
from  the  third  phase,  which  took  place  in  September,  were  used  in  this 
thesis.  A more  complete  description  of  GATE  is  available  in  Kuettner  et 
al_.  (1974). 

3 . 2 Analysis  of  Synoptic-Scale  Data 

Detailed  analysis  of  the  synoptic-scale  phenomena  in  which  the 
cells  being  modeled  occurred  is  beyond  the  scope  of  this  thesis.  It  is 
sufficient  to  categorize  the  general  synoptic  conditions  under  which  the 
cells  occurred.  Reed  ert  al_.  (1977)  developed  a useful  categorization 
scheme  by  applying  composite  analysis  techniques  to  the  eight  wave  dis- 
turbances that  crossed  western  Africa  and  the  eastern  Atlantic  during  phase 
III  and  the  preceding  interphase  period  of  GATE  (23  August  - 19  September). 
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They  divided  each  of  these  easterly  waves  into  eight  north-south  bands 
(see  Figure  3.2).  The  second  band  was  centered  on  the  region  of  maximum 
northerly  wind  component  west  of  the  trough;  the  fourth,  on  the  trough; 
the  sixth,  on  the  region  of  maximum  southerly  wind  component  east  of  the 
trough;  and  the  eighth,  on  the  following  ridge.  Bands  1,  3,  5,  and  7 occu- 
pied intermediate  positions,  The  positions  were  based  on  the  700-mb  flow 
along  the  path  of  the  wave  disturbance  center  (defined  as  the  position 
of  the  maximum  vorticity),  Reed  et_  aj^.  referred  to  these  bands  as  cate- 
gories or  phases,  but  in  this  thesis  they  will  be  referred  to  as  sectors. 
The  sectors  were  divided  into  seven  east-west  rows,  4°  of  latitude  deep 
with  the  center  of  the  disturbance  at  the  midpoint  of  the  middle  row. 

Figure  3.2  shows  how  this  process  was  accomplished  for  one  of  the  waves. 


Figure  3.  2.  Illustration  of  a GATE  wave  and  its  division  into  sectors. 

The  streamline  analysis  in  (a)  was  accomplished  on  fil- 
tered 700-mb  winds.  Orientation  of  the  features  defining 
the  sectors  is  shown  as  slanting,  dashed  lines.  The 
sectors  are  oriented  N-S  as  shown  in  (b),  centered  on  the 
points  where  the  defining  features  cross  the  path  of  the 
disturbance  center,  the  E-W  line  marked  with  a zero, 
(from  Reed,  et.  al. , 1977) 
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The  streamline  analysis  was  based  on  winds  which  had  been  band-pass  filter- 
ed to  reduce  waves  of  shorter  than  two  days  and  longer  than  six  days  to 
less  than  half-amplitude,  while  retaining  almost  fully  waves  of  three  to 
four  day  period.  The  east-west  line  marked  with  a zero  delineates  the  path 
of  the  disturbance  center.  Each  wave  was  divided  into  56  boxes  in  the 
manner  shown  and  the  observations  were  averaged  in  each  box  to  produce  the 
composite  wave. 

Reed  and  his  colleagues  found  that  the  waves  propagated  westward 
at  an  average  speed  of  8 m sec"^  in  the  zone  of  cyclonic  shear  to  the  south 
of  the  700-mb  easterly  jet.  The  mean  latitude  of  the  disturbance  paths 
was  12°N,  the  mean  wavelength  was  about  2500  km,  and  the  mean  period  was 
about  3.5  days.  The  maximum  low-level  convergence  and  upward  vertical 
motion  occurred  just  ahead  of  the  trough.  The  maximum  cloudiness  and  pre- 
cipitation was  found  just  south  of  the  disturbance  center  in  sectors  2 and 
3 near  the  position  of  maximum  convergence  and  upward  motion. 

The  paths  of  the  centers  of  the  disturbance  were  at  a fairly  con- 
stant latitude.  Thus,  the  B-scale  array  was  usually  located  in  the  row  of 
boxes  immediately  south  of  the  row  in  which  the  disturbance  paths  were 
centered.  Since  the  waves  were  usually  tilted  on  NE-SW  axes,  the  wave 
features  generally  passed  over  the  array  a few  hours  before  the  wave  sector 
indicated  that  they  did  (e,g.  the  troughs  usually  passed  over  the  array 
a few  hours  before  the  array  was  located  at  the  midpoint  of  wave  sector  4). 
This  effect  may  be  seen  in  the  wave  depicted  in  Figure  3,2. 

The  wave  sectors  discussed  in  this  section  will  be  referred  to  in 
Chapters  4 and  5 in  order  to  describe  the  general  synoptic  situations  under 
which  the  cells  being  modeled  occurred. 
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3.3  Radar  Data 

3.3.1  Source 

The  general  approach  to  the  collection  of  radar  data  during  GATE 
is  described  by  Hudlow  (1975).  The  radar  data  used  in  this  thesis  were 
collected  with  the  WR73  radar  belonging  to  the  M.I.T.  Weather  Radar  Re- 
search Project.  This  5.7-cm  radar  was  located  aboard  the  Gilliss,  one  of 
the  ships  in  the  B-scale  array  (see  Figure  3.1).  The  radar  generally 
sampled  once  every  15  minutes,  occasionally  less  often  when  significant 
convection  was  not  present.  Each  sampling  was  an  elevation  step  sequence 
in  which  360°  of  azimuth  were  scanned  at  a number  of  elevation  angles. 
This  provided  three-dimensional  data  for  constant  altitude  displays,  ver- 
tical cross-sections,  etc.  Weather  echoes  were  averaged  to  eliminate 
frequency  fluctuations,  range-normalized,  and  recorded  digitally  in  512 
range  bins  along  each  azimuth.  The  range  bins  were  set  up  to  cover  the 
256  km  range  of  the  radar  in  such  a way  as  to  provide  finer  range  resolu- 
tion near  the  radar.  This  method  gave  the  reflectivity  data  finer  reso- 
lution where  they  were  more  reliable  and  made  the  range  resolution  more 
consistent  with  the  vertical  and  azimuthal  resolutions,  A complete  des- 
cription of  the  processed  radar  data  from  the  Gilliss  may  be  found  in 
Austin  (1976). 

3.3.2  Handling 

Most  of  the  analysis  of  radar  data  in  this  study  involved  the  use 
of  three-dimensional  (3-D)  radar  maps.  These  maps  were  constructed  from 
the  polar  data  described  in  the  previous  section  by  superimposing  a three 
dimensional  rectangular  box  over  the  area  of  interest,  dividing  that  box 
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horizontally  and  vertically  into  a grid  of  smaller  boxes,  and  averaging 
the  reflectivity  data  observed  within  each  small  box.  The  result  was  a 
"stack"  of  constant-altitude  plan  position  indicator  (CAPPI)  maps,  each 
of  which  was  a horizontal  cut  through  the  region  at  a given  level  showing 
equivalent  radar  reflectivity  factor  (Z)  in  dBZ  organized  on  a cartesian 
plane.  Figure  3.3  is  an  example  of  a CAPPI  radar  map.  The  horizontal 
resolution  of  the  CAPPI' s was  1 km  for  both  the  line  considered  in  Chap- 
ter 4 and  the  cluster  considered  in  Chapter  5.  The  3-D  radar  maps  for 
the  line  were  composed  of  seven  layers,  each  2 km  deep,  so  that  the  maps 
covered  a vertical  depth  of  14  km.  The  3-D  maps  for  the  cluster  had  a 
vertical  resolution  of  1 km  and  were  composed  of  from  seven  to  ten  layers. 

Maximum-tops  radar  charts  were  used  to  compare  model  cel  1 s to  observed 
cells.  These  "max-tops"  charts  were  generated  by  determining  geometrically 
the  maximum  height  at  which  precipitation  was  detected  above  each  carte- 
sian grid  box.  This  was  accomplished  by  searching  through  polar  data  from 
progressively  higher  elevation  angles.  The  max-tops  charts  for  the  line 
had  a horizontal  resolution  of  4 km  while  the  charts  for  the  cluster  had 
a 2-km  resolution.  The  tops  in  both  sets  of  charts  were  recorded  to  the 
nearest  kilometer. 

3.3.3  Analysis 

[NOTE  ON  TERMINOLOGY:  In  this  thesis  the  terms  updraft  and  down- 
draft,  when  applied  to  the  model,  will  refer  to  the  inner  and  outer  cylin 
ders,  respectively,  of  the  model.  When  applied  to  radar  data,  the  term 
updraft  will  refer  to  the  intense  reflectivity  core  of  a cell  wherein 
active  convection  is  presumed  to  be  taking  place  and  will  be  understood 
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Figure  3.  3.  Example  of  a CAPPI  radar  map.  This  map  shows  part 
of  the  line  considered  in'  Chapter  4.  The  reflectivities, 
given  in  dBZ  for  1X1  km  grid  points,  were  determined 
by  averaging  over  the  layer  from  0-2  km.  The  3-D  map 
from  which  this  CAPPI  was  made  extended  to  a height  of 
14  km  and  was  composed  of  seven  of  these  CAPPIs. 
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to  be  analogous  to  the  inner  cylinder  of  the  model.  The  term  downdraft. 


when  applied  to  radar  data,  will  refer  to  the  region  surrounding  the 


active  convection  wherein  compensatory  vertical  motions  are  taking  place. 


and  will  be  understood  to  be  analogous  to  the  outer  cylinder  of  the  model 


These  terms  will  be  used  to  describe  regions,  not  vertical  motions,  with 


full  awareness  of  the  fact  that  descent  may  sometimes  occur  in  the  "up- 


draft" and  ascent  may  sometimes  occur  in  the  "downdraft".] 


It  was  necessary  during  the  course  of  this  study  to  inspect  an 


area  or  line  of  convective  activity  to  identify  individual  cells  and  de- 


termine the  size  of  their  updrafts.  The  precipitation  in  these  areas  and 


lines  was  nearly  always  spatially  continuous.  Thus,  it  was  difficult  to 


determine  the  boundaries  between  cells.  Since  the  nature  of  the  model 


requires  that  each  cell  be  a separate  entity,  cells  and  updrafts  had  to 


be  defined  in  such  a way  that  individual  cells  would  be  separated  from 


each  other. 


CAPPI  radar  maps,  similar  to  the  example  in  Figure  3.3,  were  used 


to  identify  cells  and  determine  their  updraft  areas.  If  one  assumes  that 


a peak  reflectin' ty  marks  the  center  of  a cell,  a reasonable  way  to  deter- 


mine the  extent  of  its  updraft  is  to  search  outward  from  the  peak  until  a 


sudden,  strong  decrease,  say  3 dB  (representing  a reduction  by  half),  in 


the  reflectivity  is  encountered.  This  point  would  mark  the  outer  edge  of 


the  uparaft.  Unfortunately,  such  a sharp  reflectivity  gradient  was  often 


absent  from  the  data,  and,  when  present,  often  delineated  unrealistically 


large  updrafts.  In  order  to  avoid  these  problems,  it  was  necessary  to 


adopt  somewhat  less  logical  but  more  broadly  applicable  rules  for  defining 


the  limits  of  an  updraft.  To  separate  cores  of  maximum  reflectivity  into 
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individual  cells,  it  was  assumed  that  the  reflectivity  between  cells  must 
at  some  point  be  less  than  or  equal  to  half  the  maximum  reflectivity. 

Thus,  a reflectivity  core  was  identified  as  a cell  if  its  point  of  maxi- 
mum reflectivity  could  not  be  connected  to  any  other  cell  without  crossing 
a point  where  the  reflectivity  was  three  or  more  decibels  less  than  its 
(own)  maximum.  This  left  some  flexibility  in  defining  the  extent  of  the 
updraft.  The  updraft  area  could  have  been  defined  as  the  area  covered  by 
data  points  containing  the  maximum  reflectivity.  However,  the  sampling 
problems  inherent  in  radar  observations  and  in  transposing  from  polar  to 
cartesian  coordinates  made  the  use  of  such  a restrictive  definition  unwise. 
Instead,  the  two  remaining  possibilities  for  the  definition  of  an  updraft 
were  used.  The  first  definition  assumed  the  updraft  area  to  be  that  area 
having  reflectivities  of  the  maximum  value  or  one  decibel  less  than  the 
maximum  value.  The  second  definition  assumed  the  updraft  area  to  be  that 
area  having  reflectivities  of  the  maximum  value  or  one  or  two  decibels 
less  than  the  maximum  value.  Both  of  these  definitions  are  consistent 
with  the  rule  used  to  separate  the  cells  in  that  neither  allows  the  up- 
drafts of  neighboring  cells  to  overlap. 

Figure  3.4  shows  how  the  process  worked.  Two  cores,  each  having 
a maximum  reflectivity  of  48  dBZ,  have  been  identified  and  are  shown  in- 
side squares.  Note  that  they  cannot  be  connected  without  passing  through 
a point  where  the  reflectivity  is  three  or  more  decibels  less  than  the 
maxima.  The  updraft  areas  based  on  the  first  definition  are  contained 
within  solid  lines.  In  the  upper  cell  part  of  this  area  is  not  contiguous 
with  the  reflectivity  core.  This  is  probably  a result  of  the  noisy  nature 
of  the  data.  Peculiarities  of  this  sort  appeared  often  in  the  radar  data 
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Figure  3.  4.  Example  of  the  use  of  a CAPPI  radar  map  to  identify 
cells  and  define  their  updraft  areas.  See  text  for 
explanation. 
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and  had  to  be  dealt  with  by  objectively  applying  specific,  subjectively- 
determined  rules.  The  areas  based  on  the  second  definition  are  contained 
within  dashed  lines.  Note  that  in  spite  of  the  fact  that  the  cells  have 
the  same  maximum  reflectivity,  their  updraft  areas  (as  defined)  differ 
considerably.  Since  the  map  has  a horizontal  resolution  of  1 km,  the  up- 
draft areas  can  be  found  by  simply  counting  the  blocks  within  the  outlined 
regions. 


3.3.4  Idealization 

A close  inspection  of  the  radar  data  in  Figure  3.4  reveals  a sub- 
stantial variability  in  cell  characteristics.  This  type  of  variability 
was  present  in  all  the  radar  data  for  the  line  and  the  cluster.  It  is  not 
clear  whether  it  was  due  to  actual  differences  in  cell  characteristics  or 
arose  from  sampling  similar  cells  at  different  times  in  their  life  his- 
tories. In  any  case,  this  variability  posed  problems  both  in  selecting 
appropriate  input  variables  for  the  model  and  in  making  comparisons  between 
observed  and  modeled  cells. 

In  order  to  determine  suitable  input  variables,  the  features  (i.e, 
the  line  and  cluster  of  cells)  were  idealized  into  geometrically  simple 
combinations  of  "representative"  cells.  The  idealized  features  had  the 
same  number  of  cells  as  the  observed  features,  but  the  representative  cells 
in  the  idealized  features  were  identical  with  each  other  and  geometrically 
compatible  with  the  model  (i.e.  they  were  composed  of  an  inner  cylinder 
which  was  the  active  core  or  updraft  region  of  the  cell,  and  an  outer 
cylinder  which  was  the  "environment"  or  area  where  compensating  motions 
y take  place.  This  idealization  process  is  described  in  section  3. 3. 4.1. 
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Since  the  model  would  be  simulating  a "representative"  cell,  it  was  also 
necessary  to  determine  representative  values  of  the  characteristics  of 
the  observed  cells  with  which  to  compare  model  cells.  The  process  for 
determining  representative  characteristics  for  comparison  is  presented  in 
section  3. 3. 4. 2. 

3. 1.4.1  Determination  of  Input  Variables 

Although  the  idealization  processes  for  the  cluster  and  the  line 
were  not  identical,  they  had  some  common  features  which  will  be  discussed 
in  this  section.  The  features  peculiar  to  the  idealization  processes  for 
the  line  and  the  cluster  will  be  discussed  in  sections  4. 1.2.1  and  5. 1.2.1, 
respectively. 

Because  it  was  necessary  to  build  some  rather  arbitrary  definitions 
into  the  idealization  process,  some  flexibility  was  also  built  into  the 
process.  As  a result  of  this  flexibility,  the  idealization  process  will 
produce  a number  of  sets  of  input  variables  from  which  a choice  of  the 
most  realistic  values  of  the  variables  may  be  made  based  on  intuitive 
judgment,  theoretical  reasoning,  and  field  observations  of  these  variables. 

The  first  step  in  the  idealization  process  was  to  determine  which 
cells  would  be  used  in  the  process.  The  method  of  identifying  cells  on  a 
CAPPI  radar  map  was  given  in  section  3.3.3.  For  the  two  CAPPI's  nearest 
the  ground,  the  cells  were  identified  and  their  maximum  reflectivity  fac- 
tors were  recorded.  Cells  were  then  said  to  "exist"  for  the  purpose  of 

the  idealization  if  they  met  certain  criteria.  The  criteria  varied  (this 

1» 

variation  provided  some  of  the  flexibility  built  into  the  process),  but  in 

\ 

all  cases  required  that  an  identified  cell  exist  in  both  of  the  lower 
layers  and  that  the  maximum  reflectivities  exceed  given  thresholds  in  one 


or  both  layers.  These  criteria  will  be  referred  to  as  the  E (for  exist- 
ence) criteria.  The  E will  be  followed  by  two  subscripts.  The  first  sub- 
script will  prescribe  the  minimum  reflectivity  required  in  one  layer  and 
the  second  will  prescribe  the  minimum  reflectivity  required  in  the  other 
layer  in  order  for  a cell  to  be  included  in  the  idealization  process. 

For  example,  E3Q  -j  will  refer  to  the  criterion  requiring  that  a cell  with 
a reflectivity  of  at  least  30  dBZ  will  exist  in  one  layer  and  some  cell 
exist  in  the  other  layer.  Similarly,  E^q  will  refer  to  the  criterion 
requiring  that  cells  with  reflectivities  of  at  least  40  dBZ  exist  in  both 
layers.  It  was  also  required  that  the  cells  in  the  two  layers  exist  in 
reasonable  geographical  proximity. 

The  method,  described  in  section  3.3.3,  for  determining  the  areas 
of  updrafts  involved  counting  the  data  points  surrounding  a core  of  maxi- 
mum reflectivity  which  had  reflectivities  within  a given  range  of  the 


maximum.  The  two  ranges  used  in  this  study  form  the  A (for  "area")  cri- 
teria. A-|  will  refer  to  an  updraft  area  determined  from  data  points 
having  the  maximum  reflectivity  or  reflectivities  one  decibel  less  than 
the  maximum.  Similarly,  A2  will  refer  to  an  updraft  area  determined  from 
data  points  having  the  maximum  reflectivity  or  reflectivities  one  or  two 
decibels  less  than  the  maximum. 

Updraft  areas  based  on  both  A-j  and  A2  criteria  were  determined 
from  the  CAPPI's  for  the  lowest  two  layers.  Although  the  updraft  area  in 
the  model  is  constant  with  height,  the  updraft  areas  determined  from  the 
radar  data  were  seldom  the  same  for  both  layers,  Thus,  methods  were  de- 
vised to  combine  the  data  from  the  two  layers  into  one  area.  The  two 


devised  methods  form  the  C (for  "combination")  criteria.  CaVg  will  refer 
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to  the  updraft  area  determined  by  averaging  the  areas  of  the  updrafts  in 
the  two  layers  while  Cmax  will  refer  to  the  updraft  area  determined  by 
using  the  maximum  area  of  the  two  layers. 

The  idealization  processes  for  the  line  and  the  cluster  employed 
different  E criteria  (i.e.  the  concept  was  the  same  but  the  subscripts 
were  different).  Both,  however,  used  the  same  A and  C criteria.  The 
downdraft  area  of  the  cells  in  the  cluster  were  determined  by  an  entirely 
different  process  than  that  used  for  the  cells  in  the  line;  these  processes 
will  be  described  in  the  following  chapters. 

3. 3. 4. 2 Determination  of  Cell  Characteristics  for  Comparison  of 
Observed  with  Model  Cells 

The  model  provides  fine  time-  and  height-resolution  information 
on  such  processes  and  properties  as  entrainment,  vertical  velocity,  con- 
densation, water  content,  forcings,  etc.  However,  little  is  known  about 
these  processes  in  real  situations  because  they  are  so  difficult  to  measure. 
The  only  property  which  was  available  from  the  model  and  which  was  observed 
in  some  useful  resolution  in  the  field  was  radar  reflectivity  (the  model 
calculates  radar  reflectivity  from  the  drop  size  distribution).  The  radar 
data  had  sufficient  spatial  resolution  to  deal  with  most  precipitating 
convective  elements.  Unfortunately,  while  its  time  resolution  (15  min) 
was  sufficient  for  studying  mesoscale  and  synoptic-scale  precipitation 
features,  it  was  too  coarse  to  allow  more  than  a single  "snapshot"  of  a 
cell  taken  at  a random  time  in  its  life.  It  was  possible,  if  the  timing 
was  right,  for  the  same  cell  to  appear  on  two  consecutive  radar  maps. 
However,  it  was  rarely  possible  to  identify  these  as  the  same  cell  because 
of  the  confusion  created  by  neighboring  cells,  Therefore,  while  the  model 
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produced  a life  history  of  the  reflectivities  within  themodeled  cells,  no 
corresponding  life  history  of  observed  cells  was  available.  The  problem, 
then,  was  to  devise  methods  by  which  fine  time-  and  height-resolution 
reflectivity  data  from  a single  cell  simulated  by  the  model  could  be  com- 
pared to  somewhat  coarser  height-resolution  reflectivity  data  from  instan- 
taneous samples  collected  by  the  radar  at  random  times  in  the  lives  of  a 
number  of  cells.  To  address  this  problem,  three  cell  characteristics  were 
chosen:  the  maximum  height,  the  maximum  reflectivity,  and  the  height  of 
the  maximum  reflectivity. 

Maximum  Height 

Iri  vigorous  convective  situations  ice  crystals  are  usually  found 
in  the  upper  portions  of  the  cells.  These  slcwly-falling  particles  may 
blow  off  as  an  anvil  or  remain  as  debris  after  the  intense  lower  portions 
of  cells  have  rained  out.  In  any  case,  they  remain  aloft  for  relatively 
long  periods  of  time,  so  that  when  a number  of  cells  occur  in  the  same 
area,  they  tend  to  form  a precipitation  canopy  over  the  area.  The  height 
of  the  canopy  is  not  as  time-dependent  as  the  features  of  the  individual 
cells  and  thus  serves  as  a good  indicator  of  the  height  which  the  more 
vigorous  cells  attained.  Not  all  cells,  of  course,  reach  the  height  of 
the  canopy,  and  occasionally  very  vigorous  cells  will  penetrate  it.  How- 
ever, it  will  be  assumed  in  this  study  that  the  height  of  the  canopy  cor- 
responds to  the  vertical  extent  of  the  representative  cells.  Maximum-tops 
radar  charts  were  used  to  determine  heights  of  cells  in  this  study. 

Maximum  Reflectivity 

During  the  cell-idealization  process,  the  maximum  reflectivity  of 


* 


each  observed  cell  in  the  lowest  two  CAPPI’s  was  determined.  This  allowed 
the  construction  of  the  distribution  of  maximum  reflectivities  in  the  up- 
drafts of  the  observed  cells.  The  model  output  allowed  the  determination 
of  the  distribution  with  time  of  the  maximum  reflectivity  in  the  model 
cell  updraft.  These  distributions  were  compared.  However,  it  is  important 
to  notice  that  the  comparison  was  made  between  the  time  spectrum  of  maxi- 
mum reflectivity  for  a single  modeled  cell  and  the  spectrum  of  instanta- 
neous maximum  reflectivity  for  a large  number  of  observed  cells.  In  order 
for  such  a comparison  to  be  meaningful,  it  must  be  assumed  that  the  ob- 
served cells  were  sampled  at  random  times  during  their  lifetimes.  Further- 
more, since  the  observed  cells  had  varying  degrees  of  representativeness, 
this  comparison  must  be  viewed  as  a somewhat  subjective  indication  of  the 
model's  ability  to  simulate  the  observed  cells. 
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Height  of  Maximum  Reflectivity 


The  approach  applied  to  comparing  maximum  reflectivities,  described 
in  the  previous  paragraph,  was  also  applied  to  comparing  the  height  of  the 
maximum  reflectivity  in  the  cell  updrafts.  The  model  output  provided  these 
heights  throughout  the  life  of  the  cell.  For  the  cluster  of  cells  consi- 
dered in  Chapter  5,  the  cell  population  and  organization  was  such  that  the 
heights  could  be  determined  by  manual  inspection  of  the  CAPPI's.  However, 
because  of  the  large  number  of  cells  and  the  nature  of  their  relationships 
with  each  other,  the  application  of  this  manual  inspection  technique  was 
not  feasible  for  the  line.  Instead,  a computer  program  was  written  which 
searched  up  each  column  in  the  3-D  radar  charts  for  the  line  and  recorded 
the  height  at  which  the  maximum  reflectivity  occurred  and  that  maximum 
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value.  The  distribution  of  height  versus  the  number  of  columns  with  maxi- 
mum reflectivities  at  that  height  was  then  found  for  reflectivities  above 
given  thresholds.  There  were  two  drawbacks  to  this  automated  approach. 
First,  except  for  the  highest  thresholds,  the  distributions  probably  in- 
cluded data  from  some  columns  which  were  not  actually  located  in  the  cell 
updrafts.  Secondly,  any  tilting  of  the  updrafts  which  occurred  might 
permit  more  than  one  maximum  height  to  be  recorded  for  a single  updraft 
and  thus  tend  to  spread  the  distribution.  These  drawbacks  increased  the 
subjectivity,  already  alluded  to  in  the  previous  paragraph,  of  this  com- 
parison scheme. 

3.4  Upper  Air  Data 

3.4.1  Source 

The  upper  air  data  used  in  this  study  is  described  in  Acheson 
(1976).  It  was  gathered  by  rawinsonde  ascents  launched  every  three  hours 
from  the  seven  ships  in  the  B-scale  array.  It  consisted  of  temperature, 
humidity,  wind,  and  height  data  every  5 mb. 

3.4.2  Handling 

Temperature  and  relative  humidity  data  were  to  be  loaded  into  the 
model  at  500  m increments  from  the  surface  to  15  km.  Thus  the  extremely 
fine  vertical  resolution  of  the  data  was  unnecessary  and  somewhat  of  a 
nuisance.  The  observed  temperature  and  relative  humidity  data  were  con- 
verted to  a form  more  consistent  with  the  model  requirements  by  taking 
centered  averages  of  the  data  from  250  m below  to  250  m aboye  the  heights 
at  which  data  would  be  loaded  into  the  model.  These  mean  temperatures  and 


- 41  - 

relative  humidities  were  then  plotted  against  a linear  height  scale  exten- 
ding from  the  surface  to  15  km. 

3.4.3  Idealization 

Although  there  were  no  soundings  taken  in  the  midst  of  the  line  or 
cluster,  there  were  numerous  soundings  available  from  their  general  vicini- 
ties. These  profiles  of  temperature  and  relative  humidity,  while  having 
major  characteristics  in  common,  varied  enough  over  space  and  time  to  make 
difficult  the  determination  of  the  actual  environment  in  which  the  convec- 
tion was  taking  place.  To  avoid  this  difficulty,  idealized  profiles, 
based  on  a number  of  soundings,  were  constructed  in  such  a way  as  to  retain 
the  general  characteristics  of  the  actual  soundings  while  ignoring  the 
minor  variations.  The  first  step  in  this  idealization  process  was  to 
identify  the  important  characteristics  --  dry  or  moist  layers,  points  where 
gradients  changed  suddenly,  etc.  Since  all  profiles  did  not  have  the  same 
characteristics,  it  was  necessary  in  sane  cases  to  identify  more  than  one 
type  of  profile.  The  properties  of  the  major  characteristics  (thickness, 
height,  maximum  or  minimum  value,  etc.)  were  then  averaged  to  produce  in- 
formation from  which  an  idealized  sounding  could  be  constructed  for  each 
type  of  profile.  The  temporal  and  spatial  distributions  of  the  idealized 
soundings  were  then  investigated  to  determine  which  would  provide  the  most 
realistic  initialization  for  the  model.  The  specifics  of  the  idealization 
process  are  given  in  the  following  chapters. 
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CHAPTER  FOUR 

SIMULATION  OF  CONVECTIVE  CELLS  IN  A LINE 

A line  of  convective  showers,  well  documented  by  radar  observa- 
tions, entered  the  GATE  B-scale  array  on  the  morning  of  5 September,  1974. 
In  this  chapter,  the  line  and  its  environment  are  described.  Complex 
variations  in  the  cells  and  the  environment  are  eliminated  through  an 
idealization  process  which  produced  "representative"  cells  in  a "repre- 
sentative" environment.  Three  sets  of  experiments  are  run  with  the  ideal- 
ized input,  and  their  results  are  compared  to  observed  data  to  determine 
how  well  the  model  simulated  the  cells. 

4.1  The  Line  and  its  Environment  --  Observed  and  Idealized 

4.1.1  Description  of  the  Line  and  its  Environment 

Marks  and  Raub  (1978),  using  radar  maps  composited  from  four  ships 
in  the  B-scale  array,  identified  three  distinct  types  of  mesoscale  features 
associated  with  the  passage  of  easterly  waves  during  phase  III  of  GATE. 

The  first  type  of  feature  was  long,  broad,  slow-moving  bands  of  deep  con- 
vection associated  with  the  wave  trough.  The  second  was  the  "squall  lines" 

— fast  moving  (11-18  m sec”^)  lines  of  deep  convection  which  moved  toward 
the  southwest  and  exhibited  a behavior  somewhat  independent  of  the  synoptic- 
scale  waves.  The  third  type  of  feature  was  long,  thin,  east-west  oriented 
lines  which  move  in  the  northerly  flow  ahead  of  the  troughs  at  speeds  of 
4-6  m sec~\  The  line  considered  in  this  chapter  is  an  example  of  this 


third  type  of  feature. 


The  origin  of  the  line  was  evident  from  satellite  pictures  from 
the  previous  day  which  showed  an  extremely  long  (>  2000  km),  narrow  string 
of  clouds  approaching  the  array  from  the  north-northwest.  At  about  2200Z 
on  4 September,  shortly  after  it  moved  within  radar  range  of  the  array, 
loosely-spaced  cells  of  moderate  intensity  developed  along  the  string  of 
clouds.  The  cells  formed  a line  about  200  km  long  and  20  km  wide.  At 
this  time,  the  array  was  in  wave  sector  2,  The  line  of  cells  expanded  in 
the  westerly  direction  and  intensified  so  that  by  the  time  the 
line  entered  the  array  at  0300Z  it  had  grown  to  be  500  km  long  and  50  km 
wide.  Figure  4.1  shows  the  line  at  0830Z  (delineated  by  arrows)  as  seen 
by  the  radar  aboard  the  Gilliss.  By  this  time,  the  wave  had  progressed 
far  enough  west  to  place  the  B-scale  array  in  wave  sector  3.  The  line 
continued  to  move  south  for  the  next  few  hours,  eventually  losing  its 
identity  in  the  broad  band  of  precipitation  associated  with  the  trough. 

The  center  of  wave  sector  4 (marking  the  trough)  passed  over  the  B-scale 
array  at  about  1800Z. 

It  appears  that  this  convective  line  was  associated  with  the  inter- 
action between  a convergence  line  embedded  in  the  northerly  flow  ahead  of 
the  trough  and  the  dynamics  of  the  trough  itself.  The  first  feature  was 
able  to  produce  only  a line  of  clouds.  However,  when  the  features  got 
close  enough  to  compound  their  effects,  the  narrow  line  of  clouds  blossomed 
into  full-scale  deep  convection. 

4.1.2  Idealization  of  the  Line  and  its  Environment 

4. 1.2.1  Idealization  of  the  Line 

Three-dimensional  radar  charts  showing  the  line  at  081 5Z  and  0830Z 
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were  constructed  as  described  in  section  3,3.2.  Cores  of  maximum  reflec- 

• 4 

tivity  were  identified  and  their  updraft  areas  were  determined  on  the 
lowest  two  CARPI 's  (0-2  and  2-4  km)  using  the  A^  and  A2  criteria  (i.e., 
the  areas  extended  IdB  and  2dB  respectively  down  from  the  reflectivity 
••  maximum  --  see  section  3. 3. 4.1).  The  cores  were  then  designated  to  repre- 

sent or  not  to  represent  cells  based  on  the  following  existence  criteria: 

Ei  1 : a core  had  to  appear  in  both  layers; 

E.0  , : a core  had  to  appear  in  both  layers  and  the  core  in 

’ at  least  one  layer  had  to  have  a maximum  reflectivity 
of  40dBZ  or  greater; 

E40  40*  a core  wlt^  a maximum  reflectivity  of  at  least  40dBZ 

’ had  to  appear  in  both  layers, 

> 

This  choice  of  existence  criteria  allowed  investigation  of  all  the  cells 
(E.|  1),  the  moderately  strong  and  strong  cells  (E^q  1 ),  or  just  the  strong- 
est cells  (E40  4(J)  appearing  in  the  line  at  the  given  times.  Each  cri- 
terion yielded  a different  number  of  cells  in  the  line.  The  updraft  area 
of  each  cell  was  then  determined  in  two  ways;  by  averaging  the  updraft 

areas  of  the  two  layers  (combination  criterion  C ) and  by  using  the 

avg 

maximum  area  of  the  two  layers  (C  ), 

max 

The  objective  of  this  process  was  to  produce  an  idealized  line 
composed  of  identical  "representative"  cells.  The  updraft  area  of  the 
representative  cells  was  determined  by  averaging  the  updraft  areas  of  the 
observed  cells  for  each  set  of  criteria.  The  existence  criteria 
determined  three  sets  of  cells  whose  areas  were  averaged.  For  each  set 
of  cells,  there  were  four  different  sets  of  areas  resulting  from  the  two 
area  criteria  and  the  two  combination  criteria.  Thus,  twelve  different 
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representative  cell  updraft  areas  were  determined.  Since  the  represen- 
tative cells  were  assumed  to  be  geometrically  compatible  with  the  model, 
the  updrafts  were  assumed  to  be  cylindrical  and  an  updraft  radius  was 
calculated  from  each  of  the  twelve  updraft  areas. 

The  last  step  in  the  idealization  process  was  to  determine  the 
radius  of  the  outer  cylinder  of  the  representative  cells.  Toward  this 
objective,  lines  were  drawn  from  cell  to  cell,  connecting  all  the  cells 
in  the  line  by  the  shortest  route  (this  was  done  for  each  E criterion). 
This  continuous,  but  crooked,  line  was  measured  to  give  the  length  of  the 
idealized  line.  It  was  then  assumed  that  the  representative  cells  were 
equally  spaced  along  the  line  and  that  the  outer  cylinders  of  the  adjacent 
cells  just  touched,  but  did  not  overlap  (thus  making  the  line  continuous 
but  precluding,  as  in  the  model,  interaction  between  cells).  The  geometry 
of  the  idealized  line  then  allowed  calculation  of  the  outer  cylinder 
radius  R from 

R = r— — v 

R 2tiPlT 

where  J l is  the  length  of  the  idealized  line  and  n is  the  cell  population. 
The  process  is  shown  schematically  in  Figure  4.2.  (a)  shows  the  line 

after  the  areas  of  the  updrafts  (shown  dark)  have  been  found.  The  cells 
have  been  connected  by  a line  of  length  Z.  (b)  shows  the  idealized  line. 
The  cells  are  equally  spaced  along  the  length  Z.  The  updrafts  are  the 


same  size,  having  areas  equal  to  the  average  of  the  areas  of  the  updrafts 
in  (a).  The  outer  cylinders  are  also  equal  in  size  and  touch  but  do  not 
overlap. 
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Table  4.1  presents  the  results  of  the  idealization  process  for  all 
combinations  of  criteria  for  the  line  at  0815Z  and  0830Z. 

It  is  interesting  to  compare  the  data  in  Table  4.1  for  the  two 
times.  The  line  is  slightly  longer  at  0830Z,  but  contains  fewer  cells. 
Furthermore,  the  updraft  radii  are  consistently  larger  (by  about  10%)  at 
0830Z.  However,  the  fact  that  the  0830Z  line  is  longer  and  has  fewer 
cells  more  than  compensates  for  the  larger  updraft  radii,  so  that  the 
resulting  a's  are  smaller  than  at  0815Z, 

As  should  be  expected  from  the  definition  of  the  criteria,  r and 
a are  smaller  for  the  A-j  cases  than  for  the  corresponding  A2  cases,  and 
are  also  smaller  for  the  C cases  than  for  the  corresponding  Cmax  cases. 

Since  E-j  y and  E^q  represent  increasing  average  cell 

strength,  it  is  possible  to  investigate  whether  the  results  of  Yau's 
sensitivity  tests  on  o are  reflected  in  the  results  of  the  idealization 
process.  Since  for  a given  A and  C criterion  the  updraft  radius  changes 
by  only  a few  percent,  the  changes  in  a as  cell  strength  increases  can  be 
compared.  This  comparison  shows  a tendency,  more  obvious  in  the  0830Zdata 
than  in  the  0815Z data,  for  o to  decrease  as  the  cells  get  stronger.  This 
is  precisely  what  the  model  predicts  will  happen,  However,  this  effect  is 
probably  a result  of  the  fact  that  the  areas  of  the  cells  eliminated  by 
the  adoption  of  more  exclusive  existence  criteria  effectively  becomes  part 
of  the  downdraft  areas  of  the  remaining  cells.  Since  the  elimination  of 
cells  does  not  appreciably  change  the  average  updraft  area,  the  increase 
in  the  downdraft  area  forces  a to  decrease.  Nevertheless,  both  the  model 
sensitivity  tests  and  the  idealization  indicate  that  stronger  cells  have  a 
larger  surrounding  area  of  downdraft. 
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Although  the  criteria  used  for  idealization  are  admittedly  arbi- 
trary, they  yielded  fairly  realistic  values  for  updraft  dimensions. 

Simpson  and  Wiggert  (1969),  using  photogrametric  techniques  and  aircraft 
penetrations,  found  updraft  radii  of  oceanic  cells  near  Puerto  Rico  rang- 
ing from  0.50  to  1.15  km.  In  later  studies  Simpson  and  VJigaert  (1971) 
used  photogrametric  techniques  alone  to  derive  updraft  radii  of  0.65  to  1 .30 
km  for  Florida  cells.  Braham  (1952)  used  data  from  aircraft  penetra- 
tions during  the  Thunderstorm  Project  to  determine  likely  radii  of  .55  to 
.75  km  for  Florida  and  Ohio  thunderstorms, 

Yau  adopted  the  hypothesis  of  Asai  and  Kasahara  (1967)  that  a 
would  take  on  a value  such  that  vertical  heat  transport  would  be  maximized. 

His  calculations  showed  that  the  value  of  o which  maximized  vertical  heat  I 

transport  varied  with  the  updraft  radius  and  the  lapse  rate,  but  generally 
fell  between  0.3  and  0,4.  The  values  of  a from  the  idealization  were  con- 
sistently smaller.  This  discrepancy  might  stem  from  improper  assumptions 
in  the  idealization  scheme.  On  the  other  hand,  it  is  possible  that  the 
cells  were  not  operating  at  peak  efficiency. 

It  appears  from  the  above  considerations  that  the  most  reasonable 
combination  of  r and  o comes  from  the  criterion  for  determining  updraft 
size. 

I 

4. 1.2. 2 Idealization  of  the  Environment 

To  determine  a representative  sounding  with  which  to  initialize 
the  model,  123  soundings  of  temperature  and  relative  humidity  for  2 and  5 
September  from  the  B-scale  ships  were  plotted  as  described  in  section  3.4.2.  j 

The  2 September  data  were  used  to  beef  up  the  sample  and  were  selected 


because  the  situation  on  that  day,  at  first  glance,  appeared  to  be  similar 


to  that  on  5 September.  The  distribution  of  the  soundings  by  day  and  by 
wave  sector  is  given  in  Table  4.2. 


L 


Table  4.2. 

Distribution  of  soundings  by 

wave 

sector  and  day. 

Wave  Sector 

2 Sep 

5 Sep 

Total 

2 

6 

10 

16 

3 

10 

16 

26 

4 

17 

23 

40 

5 

17 

23 

40 

6 

1 

0 

1 
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Relative  Humidity 

A perusal  of  the  relative  humidity  soundings  made  it  obvious  that 
there  were  two  characteristic  profiles.  In  the  first,  there  was  a moist 
layer  several  kilometers  deep  near  the  surface,  a sharp  dry  layer  a few 
kilometers  deep  above  the  surface  moist  layer,  and  a moderately  moist  layer 
above  the  dry  layer.  The  second  also  had  a moist  layer  several  kilometers 
deep  near  the  surface,  above  which  there  was  a gradual  decrease  in  rela- 
tive humidity.  The  amount  of  the  decrease  varied  so  greatly  (from  as 
little  as  10%  to  as  much  as  60%)  that  it  was  difficult  to  combine  these 
profiles  into  a single  category,  Consequently,  they  were  grouped  into  two 
types,  one  relatively  moist  aloft  and  the  other  dryer  aloft.  Thus,  three 
relative  humidity  profile  types  were  defined  as  follows: 
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Type  I:  . Moist  at  low  levels,  sharply  dryer  at  middle 
levels,  moderately  moist  at  upper  levels; 

Type  II:  Moist  at  low  levels  with  a moderate  decrease 
in  humidity  in  middle  and  upper  levels; 

Type  III:  Moist  at  low  levels  with  a strong  decrease 
in  humidity  in  middle  and  upper  levels. 

The  division  between  Types  II  and  III  soundings  was  made  at  an 
average  upper-level  relative  humidity  of  54%,  which  corresponded  to  a 
"weak  spot"  in  the  distribution  of  that  variable.  Each  sounding  was 
classified  by  type,  and  the  depth  and  average  relative  humidity  of  each 
of  its  layers  was  determined.  For  Type  I soundings,  the  height  and  value 
of  the  minimum  relative  humidity  in  the  dry  layer  was  found,  and  for  Type 
II  and  III  soundings  the  minimum  relative  humidity  in  the  upper  layer  was 
determined.  Averages  of  these  characteristics  for  each  type  of  sounding 
were  calculated.  The  results  — based  on  43  Type  I,  44  Type  II,  and  36 
Type  III  soundings  --  are  given  in  Table  4.3. 

Idealized  relative  humidity  profiles  for  each  type  were  then  con- 
structed to  reflect  the  average  characteristics.  These  profiles  are  shown 
in  Figure  4.3. 

In  order  to  determine  when  these  types  tended  to  occur,  their 
distribution  by  wave  sector  was  found.  This  distribution  is  shown  in 
Table  4.4,  in  which  (a)  is  for  2 September  and  (b)  is  for  5 September,  and 
the  numbers  in  the  parentheses  are  the  percentage  of  the  soundings  of  all 
the  given  type  which  occurred  during  the  given  sector. 
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Table  4.3.  Mean  values  of  characteristic  features  of  relative  humidity  profiles. 

Range  of  values  is  shown  in  parentheses. 

Type  Low  Layer  Low  Layer  I Mid  Layer  j Mid  Layer  Mid  Layer  I Mid  Layer  Upper  Layer  Upper  Layer 


Table  4.4.  Distribution  of  relative  humidity  sounding  types  by  wave 
sector.  Numbers  in  parentheses  are  percentage  of  soundings  of  given  type 
which  occurred  in  given  sector. 


(a) 

2 September 

(b) 

5 September 

Wave 

Sector 

Wave  ! 

Sector 

Type 

2 

3 

4 5 

6 

All 

2 

3 

4 

5 

6 

All 

I 

2 

0 

0 2 

0 

4 

18 

10 

18 

3 

0 

39 

(33) 

(0) 

(0)  (11) 

(0) 

(8) 

; (so) 

(63) 

(78) 

(13) 

(0) 

(54) 

II 

4 

8 

7 14 

1 

34 

i 

i i 

1 

4 

4 

0 

10 

(67) 

(80) 

(41)  (82) 

(100) 

(67) 

(10) 

(6) 

(17) 

(17) 

(0) 

(14) 

III 

0 

2 

10  1 

0 

13 

1 

5 

4 

70 

0 

23 

(0) 

(20) 

(59)  (6) 

(0) 

(25) 

(10) 

(31) 

(4) 

(70) 

(0) 

(32) 

All  6 10  1 


10  16  23  23 


0 


72 
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The  difference  between  the  two  days  is  immediately  obvious.  The 
Type  II  sounding  dominated  on  2 September,  and  there  were  few  Type  I 
soundings,  while  on  5 September  the  Type  I sounding  dominated  and  there 
were  few  Type  II  soundings.  The  sector  distribution  shows  that  on  2 
September,  Type  II  is  strongly  preferred  in  all  sectors  shown  except 
sector  4,  where  it  finishes  a strong  second  to  Type  III,  On  5 September 
Type  I dominates  all  sectors  except  5,  where  Type  III  is  preferred. 

Temperature 

Aside  from  slight  variations  near  the  surface  and  an  occasional 
hint  of  the  tropopause  at  15  km,  the  temperature  profiles  were  strikingly 
similar.  There  was  a temptation  to  average  a handful  of  them  and  use  that 
average  as  the  idealized  temperature  profile.  However,  the  temperature 
soundings  were  idealized  using  a method  conceptually  similar  to  that  used 
in  the  relative  humidity  idealization. 

The  obvious  cormion  characteristic  of  the  temperature  soundings  was 
a smaller  lapse  rate  in  the  lower  levels  than  in  the  upper  levels.  An 
increase  with  height  in  the  lapse  rate  occurred  between  6 and  9 km  above 
the  ground  in  each  sounding.  The  change  was  sometimes  abrupt,  but  more 
often  it  was  gradual.  Above  and  below  the  change  the  lapse  rates  were 
reasonably  constant.  Since  this  characteristic  could  be  expected  to 
affect  the  growth  of  cumulus  cells,  it  was  chosen  as  the  characteristic 
on  which  to  base  the  idealization.  The  object  of  the  idealization  was  to 
construct  a temperature-versus-height  profile  composed  of  two  straight 
lines  — obviously  a simplification,  but  more  realistic  than  assuming  a 
constant  lapse  rate  as  Yau  did  in  his  tests. 
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Toward  this  objective,  the  following  three  characteristics  were 
determined  for  109  temperature  soundings  (14  of  the  soundings  used  in 
the  relative  humidity  idealization  had  insufficient  temperature  data): 

1)  the  height  of  the  point  at  which  the  lapse  rate  changed; 

2)  the  lapse  rate  below  that  height;  and 

3)  the  overall  (0-15  km)  decrease  in  the  temperature. 

Because  the  change  in  the  lapse  rate  was  often  gradual,  taking  place  over 
two  or  three  kilometers,  the  determination  of  the  height  at  which  it 
occurred  was  somewhat  subjective.  An  attempt  was  made  to  choose  the  best 
point  at  which  to  break  the  temperature  profile  into  two  straight  lines. 
There  were  some  soundings  in  which  there  was  a sudden  decrease  in  the 
lapse  rate  near  15  km,  presumably  caused  by  the  proximity  of  the  tropo- 
pause.  For  these  soundings,  the  lapse  rate  below  this  region  was  extended 
to  15  km  in  an  attempt  to  avoid  having  the  effect  of  the  tropopause  re- 
flected in  the  lapse  rate  of  the  entire  upper  half  of  the  troposphere. 

The  average  values  of  the  three  characteristics  for  soundings 
grouped  according  to  date,  wave  sector,  and  humidity  type  are  given  in 
Tables  4. 5-4. 7.  In  Table  4,5  (a),  which  shows  the  distribution  of  the 
mean  height  of  the  change  in  the  lapse  rate  by  wave  sector,  there  is  no 
evidence  of  a consistent  pattern  of  change,  Table  4,5  (b)  shows  an  overall 
tendency  for  the  height  to  increase  from  humidity  Type  I to  Type  II  to 
Type  III.  However,  the  subjective  nature  of  the  determination  of  the 
heights  and  the  fact  that  they  were  determined  to  the  nearest  kilometer 
make  the  significance  of  such  minor  variations  questionable.  In  the  5 


Table  4.5.  Mean  heights,  in  km,  of  the  change  in  lapse  rate 
for  different  wave  sectors  (a)  and  humidity  types  (b). 

(a)  Wave  Sector 

(b)  Relative  Humidity  Type 

2 3 4 5 6 All 

I II  III  All 

2 Sep  6.8  7.2  7.7  7.7  8.0  7.5 

7.0  7.6  7.4  7.5 

5 Sep  7.9  7.9  7.5  7.8  - 7.7 

7,6  7.8  8.0  7.7 

All  7.4  7.7  7.5  7.8  8.0  7.6 

7.5  7.6  7.8  7.6 

Table  4.6.  Mean  low-level  lapse  rates,  in  °C  km  , for  different 
wave  sectors  (a)  and  relative  humidity  types  (b). 


(a)  Wave  Sector 

2 

3 

4 

5 

6 

All 

2 Sep 

5.56 

5.52 

5.56 

5.55 

5,63 

5.58 

5 Sep 

5.66 

5.61 

5.55 

5.53 

-- 

5.57 

All 

5.61 

5.68 

5.59 

5.53 

5.63 

5.57 

(b)  Relative 

Humidity  Type 

I 

II 

III 

All 

5.51 

5.54 

5.69 

5.58 

5.60 

5.47 

5.56 

5.57 

5.59 

5.53 

5.61 

5.57 

Table  4.7.  Mean  overall  (0-15  km)  temperature  decrease,  in  °C,  for 
different  wave  sectors  (a)  and  relative  humidity  types  (b). 


(a)  Wave  Sector 
2 3 


(b)  Relative  Humidity  Type 


4 

5 

6 

All 

I 

II 

III 

All 

99.9 

99.5 

100.0 

99.7 

100.0 

99.4 

100.5 

99.7 

99.2 

98.5 

-- 

99,1 

99,5 

98.5 

98.7 

99.1 

99.4 

98.9 

100,0 

99.3 

99.5 

99.2 

99.3 

99.3 

September  data  in  Table  4.6  (a),  there  is  evidence  of  a small  but  consis- 
tent decrease  in  low-level  lapse  rate  as  the  wave  approaches.  The  possible 
significance  of  this  decrease  will  be  discussed  later.  Table  4.6  (b) 
shows  no  pattern  in  the  distribution  of  low-level  lapse  rate  by  humidity 
type.  The  only  consistent  pattern  evident  in  Table  4.7  is  a small  drop 
in  the  overall  temperature  decrease  as  the  wave  approaches  on  5 September. 
This  is  consistent  with  the  decrease  in  the  low-level  lapse  rate  for  the 
same  day  seen  in  Table  4.6  (a).  The  decrease  in  low-level  lapse  rate  from 
sector  2 to  sector  5 on  5 Septemberwas  0, 13°C  km"\  Multiplying  this  by 
the  mean  height  of  the  low  layer  for  the  same  day  gives  a lapse  rate  de- 
crease of  1.0°C  per  7,7  km,  approximately  half  the  change  in  the  overall 
lapse  rate.  Therefore,  the  decrease  in  the  lapse  rate  is  spread  out  over 
the  entire  15  km.  This  result  bears  out  the  observations  of  Reed  et  al . 
(1977),  who  found  relatively  cold  air  at  850  mb  below  the  trough  and  rela- 
tively warm  air  at  300  mb  above  the  trough.  Unfortunately,  the  2 September 
data  do  not  show  the  same  effect. 

In  brief,  except  for  the  slight  decrease  in  lapse  rate  as  the  wave 
approached  on  5 September,  there  were  no  consistent  patterns  in  the  temper- 
ature profiles.  Variations  in  the  mean  characteristics  were  so  small  as  to 
be  of  questionable  significance.  It  seems  equally  realistic  to  initialize 
the  model  with  the  overall  idealization,  the  idealization  based  on  a given 
day's  data,  the  idealization  based  on  the  wave  sector,  or  the  idealization 
corresponding  to  the  appropriate  relative  humidity  profile. 
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4.2  Experiments 

Three  groups  of  experiments  were  conducted  in  an  attempt  to  simu- 
late the  cells  in  the  line.  In  the  first  group  (X1-X7)  values  of  the 
input  variables  suggested  by  the  idealization  were  used,  and  each  experi- 
ment was  a sep  rate  entity.  In  the  second  group  (X8-X12)  an  attempt  was 
made  to  humidify  the  dry  layer  in  the  Type  I relative  humidity  sounding 
by  creating  successive  cells  in  the  same  atmosphere.  Each  experiment  in 
this  group  was  initialized  using  the  final  temperature  and  relative  humi- 
dity profiles  from  the  previous  experiment.  The  third  group  of  experiments 
(X13-X15)  was  an  attempt  to  include  mesoscale  effects  in  the  development 
of  the  cells. 

Comparison  of  model  cells  with  observed  cells  is  made  after  all 
the  experiments  have  been  described,  Table  4.9  summarizes  the  values  of 
the  important  input  variables.  Some  important  characteristics  of  the 
resulting  cells  are  given  in  Table  4.10. 

4.2.1  Experiments  1-7,  Simulation  of  Cells  in  the  Line 

Experiment  1 (XI) 

The  line  of  convective  cells  being  simulated  occurred  while  the 
array  was  in  wave  sector  3,  Table  4.4  (b)  shows  that  ten  of  sixteen  sec- 
tor 3 soundings  from  5 September  were  Type  I relative  humidity  soundings. 
For  this  reason,  the  Type  I relative  humidity  sounding  was  chosen  for  this 
experiment.  The  characteristics  of  the  idealized  temperature  soundings 
were  practically  the  same  for  wave  sector  3 and  Type  I situations.  The 
chosen  temperature  profile  had  a low-level  depth  of  7.5  km,  a low-level 
lapse  rate  of  5.6°C  km"\  and  a 15-km  temperature  decrease  of  100°C  which 
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yielded  an  upper-level  lapse  rate  of  7,7°C  km~\  The  resulting  sounding, 
which  employed  the  0800Z  surface  temperature  observation  from  the  Gilliss, 
is  plotted  in  Figure  4,4. 

This  first  experiment  was  an  attempt  to  model  the  average  cell  in 
the  line.  Thus,  the  inner-cylinder  radius  and  a were  determined  from  the 
E.j  i criterion  which  used  data  from  all  cells  having  maximum  reflectivity 
cores  in  the  lowest  two  CARPI ' s . An  inner-cylinder  radius  of  1.25  and  a 
a of  0.27  were  determined  by  taking  rough  averages  of  the  values  of  these 
variables  for  the  A2  criterion  from  both  the  081 5Z  and  0830Z  idealizations 
(see  Table  4.1).  The  A2  criterion  was  used  in  order  to  produce  a large 
radius  consistent  with  the  idealization  process. 

Figure  4,5  is  a time-height  diagram  of  radar  reflectivity  factor 
(in  dBZ)  in  the  inner  cylinder  of  the  model.  The  cell  was  rather  weak  and 


attained  a height  of  only  2,5  km.  It  would  be  ludicrous  to  suggest  that 
this  simulation  resembles  the  average  cell  in  the  line,  A number  of 
reasons  for  this  failure  could  be  suggested,  but  further  experiments  will 


be  presented  before  an  attempt  is  made  to  explain  the  results  of  this 
experiment. 

Experiment  2 

It  was  expected  that  the  dry  layer  in  the  Type  I relative  humidity 
profile  would  interfere  with  the  growth  of  the  convective  cells.  However, 
in  XI,  the  precipitation  reached  a height  of  only  2,5  km,  well  below  the 
dry  layer.  Experiment  2 was  run  in  order  to  determine  whether  the  dry 
layer  played  any  role  in  that  lack  of  development.  The  Type  I relative 
humidity  profile  was  altered  to  remove  the  dry  layer  by  continuing  the  83% 


relative  humidity  characteristic  of  the  low  layer  up  to  8.5  km,  the  top 
of  the  dry  layer.  It  was  then  decreased  to  57%,  the  upper  layer  relative 
humidity. 

The  results  of  this  experiment  were,  for  all  practical  purposes, 
identical  to  the  results  of  XI.  This  suggests  that  the  dry  layer  was  not 
responsible  for  the  feebleness  of  the  cell  in  XI.  Apparently  the  83% 
relative  humidity  was  insufficient  to  support  the  condensation  necessary 
to  drive  the  updraft,  and  the  resulting  cell  was  able  to  reach  only  half 
a kilometer  above  the  saturated  velocity  perturbation. 

Experiment  3 

XI  demonstrated  that  the  model  was  unable  to  reproduce  the  deep 
convection  characteristic  of  the  line  under  the  relative  humidity  condi- 
tions associated  with  the  open  atmosphere  around  the  line.  However,  radar 
maps  (see  Figure  4.1)  suggest  that  the  Type  I relative  humidity  profile 
may  not  be  realistic.  The  individual  cells  were  embedded  in  a broad  band 
of  precipitation.  The  synoptic  situation  suggests  that  the  line  probably 
existed  in  a region  of  synoptic-scale  ascent  associated  with  the  approach- 
ing trough  and  a superimposed  region  of  mesoscale  ascent  associated  with 
the  convergence  line  along  which  the  convective  line  formed.  It  could  be 
speculated  that  this  reinforced  ascent  worked  in  concert  with  precipita- 
tion detrained  from  cells  to  produce  the  band  of  weak  precipitation  in 
which  the  thunderstorms  were  embedded.  If  this  was  the  case,  the  cells 
could  be  occurring  in  a saturated  or  nearly  saturated  environment.  X3 
was  conducted  to  investigate  the  effect  of  the  environment  being  saturated. 
All  input  variables  were  the  same  as  in  XI  except  that  the  relative  huri- 


- 67  - 


dity  was  100%  at  all  levels. 

The  result  was  the  vigorous,  deep  convection  shown  in  Figure  4.6, 
a time-height  diagram  of  reflectivity  for  the  inner  cylinder.  It  is 
interesting  to  note  that  the  cell  continued  to  grow  vertically  after  it 
had  achieved  its  maximum  reflectivity.  This  was  probably  caused  by  the 
sudden  increase  in  lapse  rate  above  7.5  km. 

Experiment  4 

The  first  three  experiments  used  values  of  r and  o which  were 
determined  from  all  the  cells  in  the  line.  In  this  and  the  next  experi- 
ment, an  attempt  was  made  to  simulate  just  the  strongest  cells.  The  inner 
radius  and  a were  found  by  averaging  the  values  for  the  E^q  criter- 

ion for  both  times.  This  included  data  from  only  those  cells  with  a 
maximum  intensity  of  at  least  40  dBZ  in  both  the  lowest  two  CAPPI's.  The 
resulting  r was  close  enough  to  the  1.25  km  used  in  the  first  three  experi- 
ments to  retain  that  value.  The  resulting  o was  0.23.  All  other  input 
variables  were  the  same  as  in  XI. 

The  result,  shown  in  Figure  4,7,  differs  only  slightly  from  the 
result  of  XI  (see  Figure  4.5).  The  cell  is  slightly  more  intense  than 
that  from  XI,  but  certainly  no  more  realistic. 

Experiment  5 

In  X3,  a respectable  cell  was  produced  by  saturating  the  atmosphere 
in  the  situation  from  XI.  In  this  experiment,  the  same  thing  is  done  for 
X4.  Thus,  a saturated  atmosphere  was  assumed  with  r = 1.25  and  a = 0.23. 
The  result,  shown  in  Figure  4.8,  was  a cell  similar  to,  but  slightly 
stronger  than  the  cell  from  X3. 


Figure  4.  7.  Time-height  diagram  of  reflectivities  for  X4. 
Contours  are  in  dBZ. 
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Experiment  6 

It  seems  likely  that  the  earlier  convective  cells  assisted  the 
presumed  mesoscale  updraft  in  humidifying  the  region  around  the  line.  If 
this  was  indeed  the  case,  the  regions  toward  the  ends  of  the  line,  where 
the  cells  were  more  widely  distributed  in  space  and  time,  would  be  expec- 
ted to  be  more  humid  than  the  surrounding  environment  (i.e.,  Type  I sound- 
ing) but  less  humid  than  the  region  in  the  center  of  the  line.  This  ex- 
periment was  run  to  see  what  type  of  cell  the  model  would  produce  in  such 
a situation.  The  relative  humidity  was  assumed  to  Le  90%  at  all  levels 
and  the  other  input  variables  were  the  same  as  in  X3. 

Figure  4.9  is  a time-height  chart  of  radar  reflectivity  for  the 
resulting  cell.  It  grew  only  half  as  tall  as  the  cell  from  X3  and  had  a 
maximum  reflectivity  of  about  36  dPZ,  Both  the  assumptions  and  the  re- 
sults of  this  experiment  represent  the  "middle  ground"  between  XI  and  X3. 

Experiment  7 

In  spite  of  the  differences  in  a,  X3  and  X5  produced  cells  of  the 
same  height.  X7  was  a deliberate  attempt  to  produce  a taller  cell  by 
changing  the  inner  cylinder  radius  to  1.50  and  a to  0.17.  Although  the 
data  show  that  neither  of  these  values  is  unreasonable,  they  would  not 
be  expected  to  occur  simultaneously  because  large  inner  radii  are  associa- 
ted with  large  o.  Thus,  this  experiment  was  more  contrived  than  realistic. 
All  variables  except  r and  a were  the  same  as  in  X3. 

The  result,  shown  in  Figure  4.10,  was  a cell  that  was  considerably 
more  vigorous  than  that  from  X3,  but  only  .5  km  taller.  These  results 
seem  to  indicate  that  even  with  a saturated  atmosphere  and  helpful  hori- 
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zontal  dimensions,  the  model  can  produce  cells  no  taller  than  about  11  km. 
It  is  possible  that  taller  ce^s  could  be  produced  by  increasing  the  lapse 
rate;  however,  because  the  vertical  temperature  profiles  were  so  consis- 
tent in  the  situation  being  modeled,  it  would  have  been  unrealistic  to 
increase  the  lapse  rate  in  order  to  try  to  produce  taller  cells. 

4.2.2  Experiments  8-12,  Humidification  of  Dry  Layer  by  Convective 
Activity 

An  analysis  of  the  moisture  budget  in  XI  showed  that  the  moisture 

Q 

perturbation  added  21.1  x 10  g of  water  to  the  inner  cylinder  while  only 

Q 

1.4  x 10  g rained  out  of  both  cylinders  combined.  Thus,  considerable 
moistening  of  the  atmosphere  took  place,  It  seemed  possible  that  if  the 
final  atmosphere  from  XI  was  used  as  the  initial  atmosphere  for  another 
cell,  further  humidification  would  take  place.  If  so,  perhaps  the  atmo- 
sphere would  be  moistened  sufficiently  to  produce  a fairly  strong  cell 
by  repeating  this  process  a number  of  times.  Experiments  8-12  tested 
this  hypothesis,  Area-weighted  means  of  the  final  temperatures  and  rela- 
tive humidities  from  XI  were  used  to  initialize  X8.  The  results  of  X8 
were  used  to  initialize  X9,  and  so  on. 

Table  4.8  gives  the  moisture  added  in  the  impulse,  the  rainfall 
from  both  cylinders  during  the  course  of  the  integration,  the  maximum 
reflectivity,  and  the  maximum  height  for  experiments  1 and  8-12,  In  X8- 
X12,  precipitation  was  still  occurring  at  the  end  of  an  hour's  integration. 
Had  the  integration  been  run  longer,  some  of  this  precipitation  would  have 
evaporated  to  add  more  moisture  to  the  atmosphere  while  the  remainder  would 
have  rained  out.  In  any  case,  it  was  lost  to  the  process  because  the  model 
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does  not  allow  for  initialization  with  ongoing  precipitation.  The  relative 
humidity  profiles  below  2 km  in  X8-X12  were  similar  enough  that  it  was 
safe  to  assume  that  approximately  the  same  amount  of  moisture  was  added 
in  the  moisture  impulses. 


Table  4.8. 

Results  of 

XI  and  X8- 

-XI 2 . 

Experiment 

Moisture 
Added 
(109  g) 

Total 
Rainfal1 
(109  g) 

Maximum 

Reflectivity 

(dBZ) 

Maximum 

Height 

(km) 

XI 

21.1 

1.4 

27.8 

2,5 

X8 

* 21-22 

4.5+ 

33.1 

4.0 

X9 

II 

5 . 6+ 

33,2 

4.5 

X10 

II 

6,5+ 

33.4 

5.0 

Xll 

II 

6.8+ 

33,3 

5.5 

X12 

U 

6,0+ 

32.5 

5.5 

Figure  4.11a  shows  relative  humidity  profiles  before  and  after  XI, 
and  after  X8  and  X9.  Figure  4.11b  shows  relative  humidity  profiles  before 
XI  and  after  X10,  Xll,  and  XI 2 . Most  of  the  moisture  from  the  first  three 
cells  goes  to  humidify  the  low  layer.  The  last  three  cells  carry  much  of 
their  moisture  into  the  dry  layer.  The  fi.nal  cell  reached  a height  of 
5.5  km. 

At  first  glance,  the  process  appears  to  be  working,  albeit  slowly. 
However,  there  were  some  discouraging  changes  taking  place  in  the  stability 
of  the  atmosphere.  Overall  ( i , e . , throughout  the  depth  of  the  model),  the 
atmosphere  was  becoming  less  stable,  a consequence  of  cooling  due  to  dry 
adiabatic  ascent  in  the  upper  regions  of  the  model,  coupled  with  heating 
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Figure  4.  11(a).  Successive  relative  humidity  profiles.  Before 
XI—,  after  XI , after  X8 , after  X9 — . 
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due  to  condensation  in  the  lower  regions.  This  increased  instability, 
however,  was  not  available  to  the  cells  because  the  lack  of  moisture 
above  6 km  cut  off  their  growth,  In  the  region  in  which  the  cells  were 
growing,  an  interesting  sequence  of  changes  in  stability  occurred.  XI 
provided  latent  heat  to  a thin  layer  near  the  ground.  X8-X11  tapped  this 
instability  along  with  the  added  moisture  to  produce  relatively  strong 
reflectivities  and  large  rain  accumulations.  However,  the  maximum  reflec- 
tivities decreased  after  XI 0 and  the  rain  accumulations  started  to  decrease 
after  XI 1.  This  seems  to  indicate  that  once  the  depth  of  the  convection 
has  more  or  less  stabilized,  the  degree  of  instability  in  the  region  in 
which  it  is  occurring  decreases  with  each  succeeding  cell.  This  comes  as 
no  surprise  — convection  exists  just  for  this  purpose.  However,  it  does 
not  encourage  acceptance  of  the  hypothesis  that  a relatively  dry  atmosphere 
may  be  humidified  by  initiating  successive  moist  convective  cells  within 
it.  It  appears  that  continuing  this  process  would  result  in  further  de- 
creasing the  instability  of  the  layer  in  which  the  cells  are  occurring. 
This,  in  turn,  would  result  in  less  vigorous  cells  whose  assaults  on  the 
dry  layer  would  become  more  feeble.  It  would  be  difficult,  short  of  ex- 
pending a great  deal  of  computer  time,  to  determine  where  this  process 
would  lead.  However,  it  is  unlikely  that  sufficient  time  would  exist  to 
carry  it  through  to  possible  fruition  without  some  interaction  with  larger- 
scale  processes  taking  place, 

4.2.3,  Experiments  13-15,  Simulation  of  a Mesoscale  Process 

The  model  formulation  does  not  allow  any  explicit  inclusion  of 
mesoscale  processes.  Some  mesoscale  processes,  however,  may  be  included 
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implicitly  in  the  choice  of  input  variables.  X3  and  X5  are  examples  of 
implicit  inclusion  of  mesoscale  effects.  In  those  experiments  the  atmo- 
sphere was  assumed  to  be  saturated,  while  any  synoptic  and  mesoscale 
motions  which  might  have  contributed  to  the  saturation  were,  by  necessity, 
ignored.  In  experiments  13-15  an  effort  is  made  to  deal  implicitly  with 
these  motions. 

An  updraft  is  weakened  in  two  ways  by  descent  in  the  outer  cylin- 
der --  the  air  that  is  entrained  carries  downward  momentum  and  has  been 
dried  by  adiabatic  descent.  However,  because  the  areas  of  outer  cylinders 
are  larger  than  areas  of  inner  cylinders,  the  compensatory  motions  in  the 
outer  cylinders  are  generally  weak,  Mesoscale  ascent  assumed  in  the  region 
of  the  line  will  work  against  the  descent  in  an  outer  cylinder  tending  to 
minimize  it  and  its  effects  on  the  updraft.  To  simulate  the  effect  of  a 
mesoscale  updraft  counteracting  the  convectively  driven  downdraft  in  an 
outer  cylinder,  a may  be  given  a very  small  value.  This  produces  an  ex- 
tremely large  outer  cylinder  and  reduces  vertical  motions  in  that  cylinder 
to  near  zero,  but  does  not  directly  affect  the  inner  cylinder.  This  ap- 
proach will  be  tested  in  XI3-X15. 

Experiment  13 

In  this  experiment,  the  input  variables  were  the  same  as  in  XI 
except  that  a was  set  to  0.01,  The  resulting  cell  is  shown  in  Figured. 12. 
In  this  experiment,  much  less  downward  momentum  was  entrained  from  the 
outer  cylinder  than  was  entrained  in  XI.  However,  the  size  of  the  outer 
cylinder  made  it  a continual  source  of  relatively  dry  air  for  entrainment 
— the  detrained  moisture  was  simply  swallowed  up  by  the  outer  cylinder. 


This  dry  air  weakened  the  cell  and  shortened  its  life. 


Experiment  14 

To  avoid  the  problem  of  the  outer  cylinder  being  an  almost  unlimi- 
ted moisture  s'nk,  the  atmosphere  was  saturated  in  XI 4.  Otherwise,  the 
input  variables  were  the  same  as  in  X13.  The  result,  shown  in  Figure  4.13, 
was  the  most  vigorous  cell  generated  thus  far.  It  had  an  unrealistically 
high  maximum  reflectivity  of  over  57  dBZ  and  reached  a height  of  12  km. 

The  dual  effects  of  a continuous  supply  of  saturated  air  and  no  downward 
momentum  to  entrain  produced  a cell  of  unrealistic  vigor. 

Experiment  15 

In  an  effort  to  "tone  down"  the  cell  generated  in  X14,  the  rela- 
tive humidity  throughout  the  atmosphere  was  lowered  to  90%  in  XI 5.  The 
resulting  cell,  shown  in  Figure  4.14,  was  surprisingly  weak.  X6,  by  com- 
parison, employed  all  the  same  input  variables  except  a = 0.27  and  was 
noticeably  more  vigorous.  The  reason  for  this  is  that  in  the  most  vigorous 
portion  of  the  updraft  where  the  turbulent  entrainment  is  greatest,  the 
updraft  in  X6  had  humidified  the  environment  to  well  above  90%  relative 
humidity,  while  in  X15  the  environmental  relative  humidity  did  not  change 
from  the  initial  90%.  Apparently,  this  moisture  effect  was  more  important 
than  the  counteracting  momentum  effect. 

It  can  be  concluded  from  this  last  series  of  experiments  that  the 
effects  of  reducing  the  downdraft  to  near  ^ero  by  making  a very  small, 
are  secondary  to  the  effects,  also  brought  about  by  making  o small,  of 


Table  4.9.  Important  input  variables  for  XI -XI 5 . 


Experiment 

# 

Inner  Radius 
(km) 

a 

Relative  Humidity 
Profile 

Temperature 

Profile 

X! 

1.25 

.27 

Type  I (see  Fig.  4.4) 

corresponding 
Type  II  relati 
humidity  profi 
(see  Fig.  4.4) 

X2 

1.25 

.27 

83%  up  to  8.5  km 

57%  above  8.5  km 

II 

X3 

1.25 

.27 

Saturated  throughout 

II 

X4 

1.25 

.23 

Type  I 

II 

X5 

1.25 

.23 

Saturated  throughout 

1 

II 

X6 

1.25 

.27 

90%  throughout 

II 

X7 

1.50 

.17 

Saturated  throughout 

II 

X8 

1.25 

.27 

Result  of  XI 
(see  Fig.  4.9a) 

Result  of  XI 

X9 

1.25 

.27 

Result  of  X8 
(see  Fig.  4.9a) 

Result  of  X8 

j 

X10 

1.25 

• 27 

Result  of  X9 
(see  Fig.  4.9a) 

Result  of  X9 

Xll 

1.25 

.27 

Result  of  X10 
(see  Fig.  4.9b) 

Result  of  X10 

X12 

1.25 

.27 

Result  of  Xll 
(see  Fig.  4.9b) 

Result  of  Xll 

X13 

1.25 

.01 

Type  I 

As  in  XI 

X14 

1.25 

.01 

Saturated  throughout 

II 

XI 5 

1.25 

.01 

90%  throughout 

II 

Table  4.10.  Characteristics  of  model  cells  from  XI -XI 5 . All  results 
are  from  inner  cylinder. 


Experiment  Max  Top  Maximum  Reflectivity  Maximum 


Cumulative  I 


# 

(km) 

dBZ 

time 

(min) 

height 

(km) 

Rainfall  Rate 
, . time 
(mm  hr"  ) !(min) 

Rainfall 

(cm) 

XI 

2.5 

27.8 

19 

2.25 

0.7 

27 

0.02 

X2 

2.5 

27.8 

19 

2.25 

0.7 

27 

0.02 

X3 

11.0 

46.9 

24 

4.25 

48.6 

34 

1.31  + 

X4 

2.5 

28.2 

18 

2.25 

0,8 

22 

0.02 

X5 

11.0 

48.6 

26 

3.75 

68.4 

32 

1.77+ 

X6 

5.5 

35.9 

23 

2.75 

6.1 

33 

0.20+ 

X7 

11.5 

53.3 

29 

3.75 

121.5 

| 34 

3.24+ 

X8 

4.0 

33.1 

21 

2.75 

2.9  ■ 

] 

28 

0.06+ 

X9 

4.5 

33.2 

21 

2.75 

3 2 ' 

1 

30 

0.08+ 

XI 0 

5.0 

33.4 

21 

2.75 

3.2  : 

30 

0.09+ 

xii 

5.5 

33.3 

21 

2.75 

3.3 

1 

30 

0.10+ 

XI 2 

5.5 

32.5 

22 

2.75 

2.8  j 

] 

29 

1 

0.09+ 

X13 

3.0 

14.2 

11 

2.25 

0.04 

i 21 

0.0006 

X14 

12.0 

57.6 

31 

3.25 

173.0  i 

35 

4.81  + 

X15 

4.0 

29.6 

17 

3.25 

1.6  | 

26 

0.03 

tp?-  g ■ 
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averaging  detrained  moisture  over  an  extremely  large  outer  cylinder. 

Since  the  purpose  of  making  ° small  was  to  simulate  the  effects  of  meso- 
scale  motion  in  counteracting  the  downdraft  in  the  outer  cylinder,  these 
results  indicate  that  the  simulation  was  not  effective. 

4.2.4.  Comparison  of  Model  Output  to  Observed  Data 

Preliminary  investigation  indicated  that  the  cells  produced  in 
experiments  3 and  5 were  the  most  likely  to  compare  favorably  with  observed 
cells.  Although  cells  from  other  experiments  will  be  mentioned  occasion- 
ally, these  two  experiments  will  bear  the  bulk  of  the  comparison,  which 
will  be  made  in  accordance  with  the  procedures  suggested  in  section  3.1 .4.2. 

Maximum  tops 

Digital  radar  maps  displaying  maximum  tops  of  echoes  to  the  near- 
est kilometer  were  constructed  from  Gilliss  radar  data  at  15-minute  inter- 
vals from  0630Z  to  1001Z.  From  these  maps,  it  was  possible  to  determine 
the  trends  in  the  height  of  the  precipitation  canopy  lying  over  the  line. 
Figure  4.15  is  an  analog  max-tops  map  for  0830Z  which  may  be  compared  with 
Figure  4.1,  the  corresponding  map  of  low-level  reflectivities.  In  Figure 
4.15  light  gray  represents  max-tops  of  6 to  9 km,  dark  gray  represents  max 
tops  of  10  to  12  km,  and  black  represents  max  tops  of  13  km.  This  map, 
which  is  reasonably  representative  of  all  the  max-tops  maps,  shows  that  in 
the  regions  of  most  active  convection,  tops  are  generally  10-12  km  with 
some  towers  to  13  km.  Tops  are  generally  lower  toward  the  ends  of  the 
line  where  cells  are  more  widely  separated  (probably  in  both  space  and  time) 
and  the  convection  is  less  vigorous.  Even  in  these  areas,  the  cells  some- 
times reached  heights  of  10-11  km,  but  apparently  their  sparsity  precluded 
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the  formation  of  a canopy.  It  appears  that  the  11  km  height  achieved  in 
X3  and  X5  was  fairly  realistic.  However,  the  model  failed  to  simulate 
those  cells  with  tops  to  13  km.  X7,  which  used  questionable  input  varia- 
bles in  an  attempt  to  produce  a taller  cell,  produced  a cell  that  reached 
only  11.5  km.  The  cell  from  XI 4 was  12  km  high,  but  was  presumed  to  be 
unrealistic  because  of  its  exceptionally  high  reflectivity.  The  failure 
of  the  realistic  cells  to  reach  13  km  might  be  attributed  to  the  loss  of 
the  latent  heat  of  fusion  in  this  version  of  the  model  which  does  not  in- 
clude the  ice  phase.  Experiments  by  Yau  (1977)  which  included  the  ice 
phase  suggested  that  its  inclusion  would  cause  model  cells  to  grow  slight- 
ly taller. 

Maximum  Reflectivity 

The  spectrum  of  observed  maximum  reflectivity  of  the  cells  in  the 
convective  line  is  plotted  in  Figure  4.16  in  the  solid  line.  This  spec- 
trum includes  all  cells  used  in  the  idealization  process  for  the  E-j  -| 
existence  criterion  for  0815Z  and  0830Z  combined  --  a total  of  74  cells. 
The  dashed  line  in  the  same  figure  is  the  time  spectrum  of  maximum  reflec- 
tivity from  X3  for  the  34  minutes  during  which  its  maximum  reflectivity 
was  within  the  range  of  values  given  on  the  graph.  The  numbers  in  the 
time  spectrum  were  multiplied  by  a constant  to  make  its  sample  size  the 
same  as  that  of  the  observed  data.  Both  spectra  were  determined  from  data 
below  4 km  averaged  over  a 2 km  depth  (i.e.  for  0-2  and  2-4  km). 

Even  if  the  model  simulated  the  representative  cell  perfectly, 
these  two  spectra  would  not  be  expected  to  be  identical.  Thus,  it  is 
difficult  to  decide  how  much  of  the  discrepancy  between  the  spectra  is 


attributable  to  variations  in  the  observed  cells  and  how  much  of  it  should 
be  blamed  on  the  model's  inability  to  simulate  the  representative  cell. 

The  differences  in  the  two  highest  reflectivity  categories  can  probably 
be  attributed  to  variations  in  the  observed  cells.  The  differences  in  the 
lower  categories  could  be  caused  by  observed  cells  being  masked  by  the 
general  precipitation.  In  any  case,  considering  the  myriad  assumptions 
and  approximations  made  to  determine  the  input  variables  for  the  experi- 
ment and  the  inherent  weaknesses  in  the  model,  the  similarity  of  the  spec- 
tra is  encouraging. 

Height  of  Maximum  Reflectivity 

The  distribution  of  observed  heights  of  reflectivity  for  the  con- 
vective line  was  found  by  the  method  given  in  section  3. 3. 4. 2.  This  dis- 
tribution included  data  from  0815Z  and  0830Z  combined  and  was  determined 
for  thresholds  of  30  and  40  dBZ.  For  comparison,  the  output  from  X3  was 
investigated  to  determine  during  what  part  of  its  lifetime  its  maximum 
reflectivity  was  greater  than  or  equal  to  the  same  thresholds.  The  dis- 
tribution of  the  height  of  the  maximum  reflectivity  during  this  portion 
of  the  life  cycle  was  then  found.  The  results  are  shown  in  Table  4.11. 

The  agreement  between  the  distributions  for  the  40  dBZ  threshold 
is  fairly  good.  The  distributions  for  the  30dBZ  do  not  agree  quite  so 
well.  A possible  reason  for  this  may  be  seen  in  Figure  4.16,  Although 
the  input  variables  for  X3  were  determined  using  cells  with  reflectivity 
maxima  as  low  as  32  dBZ,  the  contribution  by  cells  with  maxima  below 40 dBZ 
was  small.  Thus,  we  should  expect  better  agreement  for  the  comparison 
employing  the  40  dBZ  threshold. 
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Table  4.11.  Distribution  of  radar-observed  and  model  (X3)  heights  of 
maximum  reflectivity.  For  radar:  % of  columns  in  which  reflectivity 
maximum  was  at  given  height  for  given  thresholds.  For  model:  % of 
time  during  which  the  reflectivity  maximum  exceeded  the  threshold  that 
the  maximum  was  located  at  the  given  height. 


Threshold(dBZ) 

Height  (km) 

0-2 

2-4 

4-6 

6-8 

Radar 

30 

56 

34 

9 

1 

40 

49 

40 

11 

0 

Model  (X3) 

30 

68 

24 

8 

0 

40 

54 

32 

14 

0 
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Although  this  evaluation  is  subjective,  it  seems  to  indicate  that 
the  model  did  a reasonable  job  of  locating  the  height  of  the  maximum  re- 
flectivity when  a suitably  strong  threshold  was  employed.  Since  the  i"ner- 
cylinder  maximum  reflectivity,  when  it  is  stationary  or  rising,  generally 
corresponds  to  the  maximum  updraft  and,  when  descending,  usually  occurs 
just  below  the  maximum  downdraft,  it  may  be  inferred  that  the  model  is 
doing  a reasonable  job  of  locating  the  most  vigorous  vertical  motion. 

Rainfall  Rate 

Surface  rainfall  was  measured  as  the  convective  line  passed  over 
the  Gilliss  at  about  1000Z.  Five  cumulative  rainfall  measurements  allowed 
the  estimation  of  rainfall  rates  at  four  times  by  finite  differencing. 

The  observed  rainfall  rate,  estimated  from  the  rates  at  the  four  tines, 
is  plotted  in  Figure  4.17  in  a solid  line.  In  the  same  figure,  rainfall 
rates  from  X3  (dashed)  and  X5  (dotted)  were  plotted  such  that  their  maxima 
coincided  with  the  maximum  for  the  observed  rate.  The  similarity  in  the 
shapes  of  the  curves  is  striking,  especially  when  one  considers  that  the 
observed  data  were  taken  as  a cell  passed  over  the  ship  while  the  model 
data  were  calculated  at  the  base  of  the  inner  cylinder  throughout  the  life 
of  the  model.  Since  it  is  not  known  how  representative  the  cell  that 
passed  over  the  ship  was,  where  it  was  in  its  life  cycle,  or  what  portion 
of  it  passed  over  the  ship,  this  comparison  is  crude  at  best.  It  indicates 
merely  that  the  observed  and  the  model  rainfall  rates  are  in  the  same 
ballpark, 
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Figure  4.  17.  Observed  and  model  rainfall  rates.  Observed 
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ficient  to  support  deep  convection  in  the  model.  Model  cells  which  com- 
pared  favorably  with  observed  cells  were  produced  by  assuming  that  the 
atmosphere  within  the  line  was  saturated.  An  attempt  was  made  to  humidify 
the  dry  layer  in  the  Type  I relative  humidity  sounding  by  initiating 
successive  cells  in  the  same  environment,  permitting  each  of  them  to 


modify  it.  Although  the  moisture  impulses  did  add  a considerable  amount 
of  moisture  to  the  atmosphere,  the  process  was  slow  and  seemed  to  be  de- 
creasing the  available  instability.  Finally,  o was  assumed  to  be  very 


small  in  an  attmept  to  simulate  the  effects  of  mesoscale  ascent  counter- 
acting the  convectively  driven  downdraft  in  the  outer  cylinder.  The 
assumption  caused  the  outer  cylinder  to  become  an  almost  unlimited  source 
or  sink  of  moisture,  which  made  the  model  extremely  sensitive  to  humidity 
and  masked  the  desired  momentum  effects, 
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CHAPTER  FIVE 

SIMULATION  OF  CELLS  IN  A CLUSTER 

A small  cluster  of  convective  cells  developed  near  the  Gilliss 
on  the  afternoon  of  1 September  1974.  In  this  chapter,  the  cluster  and 
its  environment  are  described  and  idealized.  Two  sets  of  experiments  are 
run  using  idealized  input,  and  their  results  are  compared  to  observed  data 
to  determine  how  well  the  model  simulated  the  cells, 

5 . 1 The  Cluster  and  Its  Environment  --  Observed  arid  Idealized 

5.1.1  Description  of  the  Cluster  and  Its  Environment 
Marks  and  Raub  (1978)  found  that  there  was  little  organization  in 
the  precipitation  patterns  near  the  ridges  (wave  sectors  1,  7,  and  8) 
during  GATE  phase  III.  During  the  afternoon  of  1 September  the  B-scale 
array  was  experiencing  weak  northeasterly  flow  characteristic  of  wave 
sector  1,  and  clusters  of  convective  cells  were  scattered  randomly  about 
the  array.  The  cluster  under  consideration  in  this  chapter,  shown  circled 
in  Figure  5.1,  appeared  at  about  1430Z  30  km  south  of  the  Gilliss.  It 

persisted  for  about  three  hours,  moving  toward  the  southeast  at  about  3 

-1  2 

m sec  . At  1445Z  it  covered  an  area  of  only  30  km  but  grew  to  nearly 

2 

150  km  by  1614Z,  and  was  composed  of  from  one  to  five  identifiable  cells 
at  any  given  time  during  its  life.  Because  its  scale  was  small  compared 
to  the  spacing  of  observations,  there  was  no  evidence  of  synoptic  or  meso- 
scale  features  having  any  effect  on  the  cluster. 
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5.1.2  Idealization  of  the  Cluster  and  Its  Environment 
5. 1.2.1  Idealization  of  the  Cluster 

Three-dimensional  radar  charts  showing  the  cluster  at  15  to  30 
minute  intervals  between  1445Z  and  1715Z  were  constructed  as  described  in 
section  3.3.2.  These  charts  had  a horizontal  and  vertical  resolution  of 
1 km.  Cores  of  maximum  reflectivity  were  identified  and  their  respective 
updraft  areas  were  determined  on  the  two  lov/est  CAPPI ' s ( 0- 1 and  1-2  km) 
using  the  A^  and  A2  criteria  for  determining  updraft  areas  (see  section 
3.3.3).  The  cores  were  determined  to  represent  or  not  to  represent  cells 
based  on  five  existence  criteria.  Each  criterion  required  that  a core 
exist  in  both  layers,  and  the  following  additional  requirements  on  maximum 
reflectivities  of  the  cores  were  made; 

E.j  no  additional  requirement; 

E30  1 ; at  least  30  dBZ  in  at  least  one  layer; 

E30  3Q:  at  least  30  dBZ  in  both  layers; 

E40  3g:  at  least  40  dBZ  in  one  layer,  30  dBZ  in  the  other; 

E40  4qJ  at  least  40  dBZ  in  both  layers. 

These  criteria  are  progressively  more  exclusive,  allowing  investigation  of 

all  cells,  or  just  those  cells  above  a given  intensity.  The  updraft  area 

of  each  cell  was  then  determined  by  averaging  the  updraft  areas  in  the  two 

lowest  layers  (combination  criteria  C ) and  by  using  the  maximum  of  the 

avg 

areas  in  the  two  layers  (C  ), 

max 

The  objective  of  this  process  was  to  produce  an  idealized  cluster 


1 


composed  of  identical  "representative"  cells.  The  updraft  area  of  the 
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representative  cells  was  determined  by  averaging  the  updraft  areas  of  the 
observed  cells.  The  existence  criteria  designated  five  sets  of  cells 
whose  areas  were  averaged.  For  each  set  of  cells,  the  two  area  and  two 
combination  criteria  defined  four  sets  of  updraft  areas  to  be  averaged. 
Thus,  20  diffe  ent  representative  cell  updraft  areas  were  determined. 

For  the  line  considered  in  Chapter  4,  the  size  of  the  outer  cylin- 
der of  the  representative  cell  was  determined  from  the  geometry  of  the 
line.  This  approach  would  not  work  for  the  cluster  because  its  geometry 
was  so  ill-defined.  Instead,  the  size  of  the  outer  cylinder  was  determined 
directly  from  the  radar  data  on  the  assumption  that  the  precipitation 
surrounding  the  updrafts  was  detrained  from  the  updrafts  and  that  its 
outer  edge  marked  the  "circumference"  of  an  irregularly  shaped  outer  "cyl- 
inder". On  the  CAPPI’s  for  the  two  lowest  layers,  the  total  area  of  pre- 
cipitation echo  associated  with  the  cluster  was  divided  into  portions  such 
that  each  core  had  a separate  surrounding  area.  The  divisions  were  made 
along  the  minimum  reflectivity  between  cores.  The  area  of  the  downdraft 
of  each  cell  was  then  assumed  to  be  its  surroundingprecipitation  area.  This 
area  varied  depending  on  which  A criterion  was  used.  The  A^  criterion  produced 
a smaller  updraft  and  hence  a larger  downdraft  than  the  criterion.  The 
downdraft  areas  were  then  averaged  in  the  same  manner  as  were  the  updraft 
areas  to  produce  20  corresponding  representative  cell  downdraft  areas. 

There  is  a fundamental  difference  between  this  process  for  deter- 
mining the  downdraft  size  for  cells  in  the  cluster  and  the  process  employed 
for  cells  in  the  line.  When  more  exclusive  existence  criteria  eliminated 
cells  from  the  line,  the  method  for  determining  the  downdraft  size  had 
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the  effect  of  making  the  eliminated  cells  and  the  areas  surrounding  them 
become  part  of  the  downdraft  of  the  remaining  cells.  In  the  process  for 
the  cluster,  when  a cell  was  eliminated  by  the  adoption  of  a more  exclu- 
sive existence  criterion,  its  updraft  and  downdraft  areas  were  discarded 
and  did  not  af'ect  the  areas  assigned  to  the  remaining  cells. 

Because  no  more  than  five  cells  existed  in  the  cluster  at  any  time, 
the  data  from  ten  times  between  1445Z  and  1715Z  were  combined  to  form  a 
sample  of  25  cells  for  the  1 criterion.  The  sample  size,  of  course, 
decreased  as  E criteria  became  more  exclusive.  To  make  the  representative 
cells  geometrically  compatible  with  the  model,  the  updrafts  were  assumed 
to  be  cylindrical,  the  downdrafts  were  assumed  to  be  concentrically  cylin- 
drical, and  the  radius  of  each  was  calculated  from  their  areas.  This 
yielded  an  inner  cylinder,  radius  r and  a ratio  of  radii  a for  the  repre- 
sentative cells. 

A schematic  representation  of  the  idealization  process  is  shown  in 
Figure  5.2.  (a)  shows  the  cluster  after  the  cells  have  been  identified 

and  the  associated  precipitation  divided  into  surrounding  areas  for  each 
updraft,  (b)  shows  the  idealized  cluster  composed  of  the  same  number  of 
cells,  each  with  the  average  updraft  and  downdraft  area.  The  results  of 
the  idealization  process  are  presented  in  Table  5,1, 

The  most  interesting  feature  in  the  data  in  Table  5.1  is  the  way 
the  values  of  r and  o change  as  the  E criterion  changes  because  the  E 
criteria  define  the  strength  of  the  cells  used  in  the  idealization  process. 
The  model  predicts  that  larger  radii  lead  to  stronger  cells.  Except  for 
the  increase  in  r associated  with  the  E^q  4q  criterion,  this  effect  is  not 
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evident  from  the  data.  The  variation  in  the  strength  of  the  cells  must 
be  attributed  to  some  other  effect.  The  relation  between  the  value  of 
o and  cell  intensity  can  be  investigated  by  observing  the  dependence  of 
a on  the  E criterion.  The  model  predicts  that  for  the  same  r,  a decrease 
in  a will  result  in  a stronger  cell.  This  is  indeed  reflected  in  the 
data.  A similar  effect  was  evident  in  the  data  from  the  idealization  of 
the  line,  but  in  that  case  it  was  thought  to  be  a fortuitous  by-product 
of  the  method  used  to  determine  the  size  of  the  downdraft  areas.  In  the 
case  of  the  cluster,  on  the  other  hand,  the  effect  reflected  in  the  results 
is  more  likely  to  be  real  because  of  the  fundamental  difference,  cited 
earlier,  between  the  two  methods  used  for  determining  the  downdraft  areas. 
It  is  interesting  to  note  further  that  o shows  almost  no  change  between 
^40  30  anc*  ^40  40  cnter’’a>  This  relationship  suggests  that  the  sudden 
increase  in  r (which,  for  the  other  E criteria  had  remained  fairly  constant 
while  o changed)  caused  the  cells  to  be  stronger.  However,  because  of  the 
small  sample  size,  this  must  be  regarded  as  speculation. 

The  inner cyl inders  for  the  cluster  were  smaller  than  those  for  the 
line,  but  did  fall  within  the  range  of  values  of  updraft  radii  observed 
by  Simpson  and  Wiggert  (1969,  1971),  The  values  of  a for  the  cells  in 
the  cluster  were  within  the  range  of  values  found  for  the  cells  in  the 
line. 

5, 1,2.2  Idealization  of  the  Environment 

To  determine  a representative  sounding  with  which  to  initialize 


the  model,  27  soundings  from  the  B-scale  ships  were  plotted  as  described 
in  section  3.4.2.  The  soundings,  taken  every  three  hours,  spanned  the 
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period  from  0600Z  to  1800Z  on  1 September.  The  first  few  hours  of  the 
period  were  from  sector  8 of  the  previous  wave,  but  the  bulk  of  the  data 
were  from  sector  1 of  the  approaching  wave.  Because  the  soundings  were 
not  distributed  over  a number  of  wave  sectors,  no  attempt  was  made  to 
categorize  them  by  sector.  Instead,  they  were  categorized  by  time  to 
allow  investigation  of  temporal  variability  and  construction  of  the  best 
possible  idealized  sounding  for  the  time  when  the  cluster  occurred. 

Relative  Humidity 

All  27  relative  humidity  soundings  were  moist  at  low  levels  with 
a strong  decrease  in  humidity  in  the  middle  and  upper  levels  and  could 
have  been  classified  as  Type  III  soundings.  Within  the  Type  III  regime, 
however,  the  soundings  had  strikingly  similar  characteristics.  Each  was 
composed  of  six  layers,  alternately  moist  and  dry,  which  covered  approx- 
imately the  lowest  10  km.  Above  the  sixth  layer,  the  relative  humidity 
was  higher  than  in  the  sixth  layer  and  reasonably  constant. 

To  develop  idealized  soundings,  each  layer  was  characterized  by 
its  height  and  relative  humidity  and  averages  of  these  quantities  were 
calculated  for  each  time  and  for  all  the  data  combined.  The  moist  layer 
immediately  above  the  surface  was  defined  by  its  surface  relative  humidity, 
the  location  and  value  of  its  maximum  relative  humidity,  and  the  height 
and  relative  humidity  at  its  top.  The  next  five  layers  were  defined  by 
the  location  and  value  of  the  minimum  (for  the  dry  layers)  or  maximum  (for 
the  moist  layers)  relative  humidity.  The  region  above  the  sixth  layer 
was  defined  by  the  height  of  its  base  and  its  mean  relative  humidity. 

The  results  of  this  process  are  shown  in  Table  5.2.  Figure  5.3, 
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Table  5.2.  Mean  values  of  characteristic  features  of  the  idealized  relative  humidity  profiles. 
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a plot  of  the  idealized  profile  for  all  the  data,  shows  the  general  shape 
of  the  profiles.  Table  5.2  shows  that  the  mean  heights  of  the  features 
were  consistent  from  time  to  time,  seldom  varying  by  more  than  a half 
kilometer  and  never  varying  by  more  than  one  kilometer.  The  minimum 
relative  humidity  in  each  dry  layer  was  less  than  that  in  the  dry  layer 
below  it,  and  the  same  was  true  for  the  maximum  relative  humidities  in 
the  moist  layers.  As  the  day  proceeded,  the  vertical  variations  in  humi- 
dity became  less  sharp  the  dry  layer  minima  increased  and  the  moist 
layer  maxima  decreased.  The  top  layer  showed  a significant  increase  in 
relative  humidity  throughout  the  12-hour  period. 

Temperature 

The  temperature  profiles  associated  with  the  cluster  were  similar 
to  each  other  and  displayed  the  same  characteristics  as  those  associated 
with  the  line:  a stronger  lapse  rate  aloft  than  near  the  surface  and  a 
point  between  6 and  9 km  where  the  change  in  the  lapse  rate  occurred. 
Therefore,  idealization  of  temperature  profiles  for  the  cluster  was  accom- 
plished in  the  same  manner  as  for  the  line  (see  section  4. 1.2. 2)  except 
that  the  data  were  categorized  by  time  instead  of  by  wave  sector  and 
relative  humidity  type.  The  results,  shown  in  Table  5.3,  indicate  that 
there  was  little  variation  in  the  temperature  profile  during  the  course 
of  the  day. 

5.2  Experiments 

Two  groups  of  experiments  were  conducted  in  an  attempt  to  simulate 


Table  5.3.  Mean  values  of  characteristic  features 
of  temperature  profiles  for  the  cluster. 

Time  Sample  size  Height  of  change  Low  level  lapse  Overall  (0-15  km) 
(Z)  in  lapse  rate  (km)  rate  (°C/km)  temp  decrease  (°C) 


99.40 

98.83 

98.75 

98.20 


99.00 


06-18  27 


98.84 
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the  cells  in  the  cluster.  The  general  approach  to  these  experiments  was 
the  same  as  the  approach  used  in  the  experiments  for  the  line,  and  most 
of  the  experiments  in  this  chapter  correspond  to  similar  experiments  in 
Chapter  4.  The  first  group  of  experiments  (X16-X19)  employed  values  of 
the  input  variables  suggested  by  the  idealization.  In  the  second  group 
(X20-X23)  an  attempt  was  made  to  humidify  a dry  layer  by  initiating  suc- 
cessive cells  in  the  same  atmosphere.  Each  experiment  in  this  group  was 
initialized  using  the  final  temperature  and  relative  humidity  profiles 
from  the  previous  experiment.  Table  5,5  summarizes  the  values  of  the 
important  input  variables  in  the  experiments.  Some  important  character- 
istics of  the  resulting  cells  are  given  in  Table  5.6. 

5.2.1  Experiments  16-19,  Simulation  of  Cells  in  the  Cluster 

Experiment  16  (XI 6) 

The  temperature  and  relative  humidity  soundings  with  which  XI 6 
was  initialized  were  the  idealized  profiles  for  1500Z,  the  only  time  at 
which  soundings  were  taken  during  the  life  of  the  cluster.  The  tempera- 
ture profile  employed  the  surface  air  temperature  from  the  1500Z  Gilliss 
observation.  The  sounding  is  plotted  in  Figure  5,4. 

This  experiment  was  an  attempt  to  model  the  average  cell  in  the 
cluster.  Thus,  the  inner  cylinder  radius  and  a we re  determined  from  the 
E^  j criterion  which  used  data  from  all  cells  having  reflectivity  cores 
in  the  lowest  two  CAPPI's.  A relatively  large  radius  was  chosen  by  using 
the  A,  C w criteria  which  yielded  a radius  of  0.8  km  ani  a o of  0,30 
(see  Table  5.1). 

Figure  5.5  is  a time-height  diagram  of  the  resulting  radar  reflec 
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tivity  factor  (in  dBZ)  for  the  inner  cylinder  of  the  model.  The  cell 
was  not  very  strong  and  attained  a height  of  only  3 km.  However,  it 
should  not  be  written  off  as  unrealistic  as  was  the  cell  from  the  corres- 
ponding experiment  (XI)  for  the  line  because  the  cells  in  the  cluster 
appeared  to  be  quite  different,  in  some  respects,  from  the  cells  in  the 
line. 


Experiment  17 

In  this  experiment  the  initial  relative  humidity  throughout  the 
atmosphere  was  set  to  100%  under  the  assumption  that  the  environment  had 
been  saturated  by  previous  convective  activity  and  mesoscale  ascent 
(when  the  assumption  was  applied  in  Chapter  4,  the  most  realistic  cells 
were  produced).  All  other  input  variables  were  the  same  as  in  X16. 

The  resulting  cell,  shown  in  Figure  5.6,  grew  to  a height  of  10  km 
and  had  a maximum  reflectivity  of  nearly  45  dBZ. 

Experiment  18 

The  mesoscale  effects  associated  with  the  cluster  were  probably 
weak  and  short-lived  compared  with  those  associated  with  the  line.  Further- 
more, the  convective  activity  had  not  been  occurring  for  nearly  so  long  a 
time.  It  is  possible  that  although  these  effects  worked  to  humidify  the 
cluster  environment,  they  were  £00  weak  and  short-lived  to  saturate  it. 

This  experiment  investigated  the  type  of  cell  which  would  grow  in  a humidi- 
fied, but  not  saturated,  environment.  The  relative  humidity  was  assumed 
to  be  90%  throughout  the  depth  of  the  model  and  all  other  input  variables 
were  the  same  as  in  XI 7. 

Figure  5.7  shows  the  resulting  cell.  It  is  essentially  a scaled- 


- 113  - 


r) 


down  version  of  the  cell  from  XI 7 . It  achieved  a height  of  6.5  km  and  a 
maximum  reflectivity  of  36.4  dBZ. 

I L_1 

Experiment  19 

Experiments  16-18  employed  values  of  r and  a based  on  all  the  cells 

I’ 

in  the  cluster.  In  this  experiment,  those  variables  were  chosen  from  the 

^40  40  ^2  Snax  cr"’teria  (r  = °-90,  cr  = 0.21)  to  simulate  the  strongest 
cells.  The  relative  humidity  was  assumed  to  be  90%  throughout  the  atmo- 
sphere and  the  temperature  profile  was  the  same  as  that  used  in  the  pre- 
vious experiments.  Figure  5.8  is  the  time-height  diagram  of  reflectivity 
for  the  resulting  cell.  The  cell  grew  no  taller  than  that  from  X17,  but 
achieved  a maximum  reflectivity  of  nearly  40  dBZ. 

5.2.2  Experiments  20-23,  Humidification  of  Dry  Layer  by  Convective 
Activity 

In  Chapter  4,  a series  of  experiments  was  described  that  tested 
the  hypothesis  that  successive  cells  initiated  in  the  same  environment 
would  progressively  humidify  the  environment  and,  in  turn,  would  result 
in  successively  stronger  cells.  Although  the  hypothesis  was  not  disproved, 
stability  and  time  considerations  made  it  questionable.  The  idealized 
relative  humidity  profile  for  the  cluster  was  similar  to  that  for  the  line 
in  that  both  contained  dry  layers  sufficiently  close  to  the  surface  to 
prematurely  disrupt  the  convective  process.  However,  the  profile  for  the 
cluster  was  dissimilar  enough  to  warrant  further  evaluation  of  the  hypo- 
thesis. The  process  used  in  this  section  was  the  same  as  that  used  in 
Chapter  4.  Briefly,  area-weighted  averages  of  the  final  temperature  and 
relative  humidity  profiles  of  both  cylinders  from  XI 6 were  used  to  initial- 
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diagram  of  reflectivities  for  X 19. 
: in  dBZ. 
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ize  X20,  the  averages  from  x20  were  used  to  initialize  X21,  and  so  on. 

Table  5.4  gives  the  moisture  added  in  the  impulse,  the  rainfall 
from  both  cylinders  during  the  course  of  the  integration,  the  maximum 
reflectivity,  and  the  maximum  height  for  experiments  16  and  20-23.  It 
is  obvious  from  a comparison  of  the  moisture  added  and  the  total  rainfall 
that  the  environment  was  being  humidified.  Figure  5,9  shows  the  succes- 
sive relative  humidity  profiles. 


Table  5.4. 

Results 

of  XI 6 and 

X20-X23. 

Experiment 

Moi sture 

(109  g) 

Total 

Rainfall 

(108  g) 

Maximum 

Reflectivity 

(dBZ) 

Maximum 

Height 

(km) 

XI 6 

7.8 

3.0 

26.7 

3.0 

X20 

9.7 

4.4 

29.3 

3.5 

X21 

10.9 

8.1 

32.4 

4.0 

X22 

11.7 

12,7 

33.3 

4.5 

X23 

12.2 

18.5 

35.3 

5.0 

The  process  showed  more  promise  for  the  cluster  than  it  did 
for  the  line.  In  the  line  situation  the  dry  layer  was  extremely  sharp. 

The  added  moisture  was  insufficient  to  humidify  it  very  rapidly  with  the 
result  that  the  capped  convection  occurred  continually  within  the  same 
layer  and  began  destroying  its  instability.  The  successive  cells  then 
became  weaker,  and  it  seemed  unlikely  that  the  process  would  lead  to  cells 
of  the  same  strength  and  vertical  extent  as  those  actually  observed  in 
the  line.  In  the  idealized  relative  humidity  profile  for  the  cluster,  the 
decrease  in  humidity  in  the  dry  layers  was  more  gradual.  The  lowest  dry 
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layer  was  quickly  humidified  by  the  first  three  cells  and  the  last  two 
cells  raised  the  second  dry  layer  a half  kilometer.  The  cells  were  grow- 
ing progressively  stronger  and  taller.  It  is  unlikely,  however,  that  such 
progress  would  continue  to  be  made  because  the  process  had  produced  a 
sharp  dry  layer  similar  to  that  encountered  in  the  line  situation  (compare 
X10,  Figure  4.11  (b)  to  X23,  Figure  5,9  (b)).  Thus,  one  might  expect  that 
the  instability  in  the  layer  in  which  the  convection  was  occurring  would 
be  weakened  by  the  next  few  cells,  and  the  process  would  suffer  from  the 
same  problems  which  arose  in  the  line  situation.  Why,  then,  did  this 
process  show  more  promise  for  the  cluster  than  it  showed  for  the  line? 
There  are  two  reasons.  First,  the  relative  humidity  profile  was  such  that 
the  successive  cells  could  grow  steadily  in  height  resultin'"  in  increasing 
instability  from  one  cell  to  the  next.  This  allowed  successive  cells  to 
become  more  intense.  Secondly,  the  cells  observed  in  the  cluster  were, 
on  the  average,  not  nearly  as  tall  or  intense  as  the  ones  observed  in  the 
line.  Thus,  by  the  time  the  effects  of  the  evolved  sharp  dry  layer  could 
be  felt,  the  model  cells  had  become  fairly  realistic, 

One  major  drawback  of  the  hypothesized  process  as  a means  for 
humidifying  dry  layers  is  that  it  proceeds  too  slowly,  as  seen  in  both  the 
line  and  cluster  situations,  It  is  possible  that  a cell  could  be  initia- 
ted before  the  previous  cell  had  completely  dissipated,  allowing  the 
process  to  run  faster  (the  formulation  of  the  model  does  not  allow  for 
the  testing  of  this  possibility).  The  stability  problems,  however,  would 
remain. 


2 
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Table  5.5.  Important  input  variables  for  X16-X23. 


Experiment  Inner  Radius  ! 0 Relative  Humidity  Temperature 

# (km)  I Profile  Profile 


15Z  (see  Fig.  5.4);  157.  (see  Fig.  5.4 


100%  throughout 


90%  throughout 


Result  of  XI 6 
(see  Fig.  5.9(a)) 


Result  of  XI 6 


Result  of  X20 
(see  Fig.  5.9(a)) 

Result  of  X21 
(see  Fig.  5.9(a)) 

! Result  of  X22 
j (see  Fig.  5.9(b)) 


Result  of  X20 


Result  of  X21 


Result  of  X22 


5.2.3  Comparison  of  Model  Output  to  Observed  Data 

The  general  approach  to  this  comparison  will  be  that  described  in 
section  3. 3. 4. 2.  However,  an  approach  geared  more  toward  individual  cells 
than  the  broad  statistical  approach  used  in  Chapter  4 will  be  applied  to 
the  cluster  because  the  cells  observed  in  the  cluster  were  less  numerous 
and  easier  to  identify  in  all  the  CAPPI  layers.  Furthermore,  these  cells 
showed  considerably  greater  variability  than  those  in  the  line,  thus 
making  a more  individualized  approach  to  their  comparison  more  meaningful. 

Maximum  Tops 

The  maximum  tops  of  the  cells  in  the  cluster  were  not  nearly  as 
consistent  as  they  were  for  the  line,  and  there  was  little  evidence  of  a 
precipitation  canopy.  Early  in  the  life  of  the  cluster,  when  the  cells 
were  weak,  they  generally  reached  2-3  km.  Later,  when  the  cells  had 
reached  their  maximum  strength,  they  extended  to  6-7  km.  Initially,  the 
observed  tops  were  consistent  with  the  results  of  X16.  This  suggests  that 
the  dry  layers  evident  in  the  large-scale  environment  existed  in  the  imme- 
diate vicinity  of  the  cluster  as  it  began  to  grow;  i.e.,  there  was  probably 
no  significant  modification  of  the  humidity  profile  prior  to  the  onset  of 
cell  growth.  The  tops  associated  with  the  strongest  observed  cells  were 
consistent  with  the  results  of  XI 8 and  XI 9,  which  suggests  that  the  cell 
environment  within  the  cluster  had  been  humidified  to  some  extent.  The  10 
km  top  produced  in  XI 7 was  well  in  excess  of  observed  tops,  which  suggests 
that  the  assumed  100%  relative  humidity  was  too  high  or  was  extended  too 
high  in  the  atmosphere. 


Maximum  Reflectivity 

The  spectrum  of  observed  maximum  reflectivity  of  cells  in  the 
cluster  is  plotted  in  Figure  5.10  in  the  solid  line.  This  spectrum  includes 
the  25  cells  used  in  the  idealization  process  of  the  ^ existence  criterion. 
The  dashed  line  in  the  same  figure  is  the  time  spectrum  of  maximum  reflec- 
tivities from  the  cel  1 produced  in  XI 8 for  the  45  minutes  during  which  its 
maximum  reflectivity  was  within  the  range  of  values  given  on  the  graph. 

The  numbers  in  the  time  spectrum  were  multiplied  by  a constant  to  make  its 
sample  size  the  same  as  that  of  the  observed  data.  In  the  discussion  in 
Chapter  4 it  was  noted  that  if  it  is  assumed  that  the  cells  were  observed 
at  random  times  in  their  lives,  the  differences  between  these  two  types 
of  spectra  can  be  attributed  to  variations  in  the  characteristics  of  the 
observed  cells  or  to  the  inability  of  the  model  to  simulate  the  "repre- 
sentative" cell.  The  reflectivity  of  the  model  cell  did  not  exceed  35  dBZ, 
whereas  over  one-third  of  the  observed  cells  had  reflectivities  greater 
than  that  value.  This  discrepancy  is  probably  due  to  the  large  variation 
in  the  characteristics  of  the  observed  cells  which,  in  turn,  may  reflect 
variations  in  their  environments.  One  model  cell  simply  could  not  be  ex- 
pected to  embody  the  characteristics  of  such  a wide  variety  of  observed 
cells.  To  make  the  comparison  somewhat  more  realistic,  the  observed  cells 
having  reflectivities  greater  than  35  dBZ  were  removed  from  the  sample. 

The  time  spectrum,  multiplied  by  a different  constant  to  make  its  sample 
size  agree  with  the  new  observed  sample  size,  is  shown  in  the  dotted  line. 
It  agrees  much  better  with  the  observed  data  than  the  first  time  spectrum. 
This  result,  rather  than  indicating  a limitation  of  the  model,  points  to  a 
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weakness  in  the  idealization  process:  it  is  not  reasonable  to  expect  a 
"representative"  cell  to  be  very  representative  of  a number  of  cells  with 
highly  variable  characteristics. 

A more  subjective,  individualized  approach  to  comparing  the  maxi- 
mum reflectivities  of  the  observed  and  model  cells  suggests  that  the  model 
cells  were  not  unrealistic.  The  maximum  reflectivity  of  the  cell  produced 
in  XI 6 is  similar  to  maximum  reflectivities  of  cells  observed  early  in 
the  life  of  the  cluster.  The  cells  generated  in  XI 8 and  X19  resembled 
cells  observed  later  in  maximum  reflectivity  in  the  same  way  that  they 
resembled  those  cells  in  maximum  height,  This  strengthens  the  suggestion 
made  in  the  previous  section  that  the  cell  environment  within  the  cluster 
was  humidified  during  the  life  of  the  cluster.  Experiment  17  produced  a 
maximum  reflectivity  of  about  45  dBZ,  representative  of  the  strongest  cells 
in  the  cluster.  Its  maximum  top  of  10  km  was  much  higher  than  observed 
maximum  tops;  however,  this  feature  could  probably  be  controlled  without 
significantly  affecting  the  maximum  reflectivity  by  decreasing  the  rela- 
tive humidity  above  7 km. 

Height  of  Maximum  Reflectivity 

The  simplicity  of  the  radar  data  for  the  cluster  allowed  each  cell 
to  be  identified  in  each  CAPRI  layer.  Therefore,  it  was  not  necessary  to 
resort  to  the  statistical  method  described  in  section  3, 3. 4. 2 to  determine 
the  distribution  of  heights  of  the  reflectivity  maxima.  That  distribution 
for  the  cell  from  X18  was  found  explicitly  by  first  determining  during 
what  part  of  its  life  its  maximum  reflectivity  exceeded  18  dBZ,  the  lowest 
maximum  reflectivity  of  any  cell  in  the  cluster.  Then  the  distribution  of 


- 126  - 


of  the  height  of  the  model  maximum  reflectivity  during  this  portion  of 
the  cell  life  was  found.  The  observed  and  model  distributions  are  shown 
in  Table  5,7.  No  cell  was  observed  to  have  its  maximum  reflectivity  above 
2 km,  but  the  model  cell  had  its  maximum  as  high  as  3-4  km.  This  could  be 
attributed  in  part  to  the  many  observed  cells  which  were  not  very  tall  and 
thus  would  tend  to  shift  the  distribution  of  maximum  reflectivity  heights 
toward  lower  values.  However,  there  were  enough  strong  cells  that  some 
higher  values  should  have  been  observed  if  they  had  existed.  Therefore, 
there  seems  to  be  a discrepancy  between  the  model  and  observed  cells  in 
this  respect. 


Table  5.7.  Distribution  of  radar-observed  vs, 
model  (X18)  heights  of  maximum  reflectivity. 
For  radar:  % of  cells  in  which  reflectivity 
maximum  was  located  at  the  given  height.  For 
model:  % of  time  during  which  the  reflectivity 
maximum  exceeded  18  dBZ  that  the  maximum  was 
located  at  the  given  height. 


Height  (km) 

0-1 

1-2 

2-3 

3-4 

Radar 

75 

28 

0 

0 

Model 

(X18) 

62 

7 

17 

14 
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5.3  Summary 

The  cells  in  a cluster  which  occurred  on  1 September  in  the  GATE 
B-scale  array  have  been  modeled  using  input  variables  determined  by  ideal- 
izing the  cluster  and  its  environment.  The  moisture  in  the  open  atmo- 
sphere around  the  cluster  was  found  to  be  insufficient  to  support  in  the 
model  the  vigorous  convection  associated  with  the  strongest  cells  observed 
in  the  cluster.  However,  the  open-atmosphere  moisture  profile  did  produce 
model  cells  similar  to  those  in  the  cluster  at  the  beginning  of  its  life 
when  the  convection  was  not  as  vigorous.  The  model  produced  cells  which 
compared  favorably  with  the  more  vigorous  observed  cells  when  the  environ- 
ment of  the  cluster  was  assumed  to  have  been  humidified  to  90%  relative 
humidity.  An  attempt  was  made  to  humidify  the  dry  portions  of  the  atmo- 
sphere near  the  ground  by  initiating  successive  cells  in  the  same  environ- 
ment and  allowing  each  to  add  moisture  to  it.  Although  this  process  was 
more  successful  for  the  cluster  than  it  was  for  the  line,  time  and  stabili- 
ty considerations  made  it  seem  unlikely  that  the  process  could,  working 
by  itself,  humidify  a sufficient  depth  of  the  atmosphere  rapidly  enough 
to  allow  the  model  to  produce  cells  having  the  strength  and  height  of  the 
more  vigorous  observed  cells,  It  was  recognized  that  the  large  variability 
of  the  cell  characteristics  in  the  cluster  made  the  "representative"  cell 
determined  by  the  idealization  process  actually  representative  of  only 
a small  percentage  of  the  cells  in  that  cluster. 
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CHAPTER  SIX 
ASSESSMENT  OF  RESULTS 

The  model  has  been  tested  against  cells  observed  in  two  situations. 
In  the  first,  the  cells  occurred  in  a line  which  was  well -developed  and 
long-lived.  This  convective  line  occurred  during  a synoptically  active 
period  when  an  easterly  wave  was  approaching  the  GATE  B-scale  array.  The 
second  situation,  on  the  other  hand,  was  relatively  undisturbed;  a synop- 
tic-scale ridge  was  in  the  vicinity  of  the  array.  In  this  situation,  the 
cells  occurred  in  a cluster  which  was  short-lived  and  poorly  developed. 

In  order  to  initialize  the  model,  it  was  necessary  to  specify  the 
horizontal  dimensions  of  the  cell  being  modeled  and  the  temperature  and 
relative  humidity  profiles  of  the  atmosphere  in  which  it  was  occurring. 
Because  of  the  variations  observed  in  cell  characteristics  and  the  fact 
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that  the  observed  cells  were  geometrically  more  complex  than  the  model 
formulation,  it  was  necessary  to  idealize  the  numerous  cells  in  the  fea- 
tures to  determine  the  dimensions  of  a "representative"  cell.  The  ideali- 
zation process  was  somewhat  arbitrary,  but  produced  updraft  dimensions 
comparable  with  observed  values  and  downdraft  (i.e.  outer  cylinder)  dimen- 

I 1 

sions  which  seemed  reasonable  based  on  theoretical  reasoning  and  intuitive 
judgment. 

There  were  no  soundings  available  from  within  the  mesoscale  fea- 
tures containing  the  cells.  Soundings  from  the  open  atmosphere  around 
those  features  had  major  characteristics  in  conmon,  but  were  not  identical. 

They  were  idealized  to  retain  the  major  characteristics  while  smoothing 
out  the  minor  variations.  The  relative  humidity  soundings  fell  into  three 
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idealized  types,  the  distribution  of  which  varied  with  the  situation. 

On  the  other  hand,  there  was  only  one  type  of  temperature  sounding,  and 
the  variations  in  its  characteristics  were  insignificantly  small  and  not 
systematic.  The  idealization  process  was  important  to  the  determination 
of  a representative  relative  humidity  profile  because  it  retained  the 
important  characteristics  which  probably  would  have  been  lost  or  signi- 
ficantly altered  in  a mean  relative  humidity  sounding.  However,  because 
of  the  lack  of  variation  in  the  temperature  profiles,  a mean  temperature 
sounding  would  probably  have  served  as  well  as  the  idealized  soundings 
used  in  the  experiments. 

In  the  case  of  the  line,  the  model  cells  initialized  with  the 
relative,  humidity  profile,  idealized  from  open-atmosphere  soundings  around 
the  line,  were  unreal istically  weak.  However,  the  model  generated  cells 
which  were  quite  similar  to  the  observed  cells  when  some  moisture  was 
added  to  the  atmosphere.  The  addition  of  moisture  was  justified  on  the 
assumption  that  the  cells  were  growing  in  a mesoscale  region  which  was  more 
moist  than  the  open  atmosphere.  The  results  of  further  experiments  sugges- 
ted that  the  additional  humidity  required  for  realistic  model  cells  could 
not  be  added  solely  by  the  action  of  previous  convective  cells.  Apparent- 
ly, some  mesoscale  convergence  and  ascent  was  also  needed. 

In  the  case  of  the  cluster,  the  model  cell  based  on  the  open-atmo- 
sphere relative  humidity  did  not  compare  favorably  with  the  representative 
cell,  but  was  similar  to  the  cells  which  occurred  early  in  the  life  of  the 
cluster.  As  in  the  case  of  the  line,  it  was  necessary  to  add  moisture  to 
produce  cells  similar  to  the  representative  cell  in  the  cluster.  Although 
further  calculations  showed  that  previous  convective  activity  could  con- 
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tribute  significantly  to  the  humidification  of  the  atmosphere,  this  con- 
tribution was  rather  small  and  would  require  an  unrealistically  long  time 
to  humidify  the  environment  to  the  point  where  representative  cells  could 
be  produced.  Once  again,  this  result  suggested  that  some  mesoscale  ascent 
was  necessary. 

The  techniques  for  deducing  representative  cell  dimensions  from 
radar  data  appear  to  be  satisfactory.  However,  in  this  study  the  process 
was  accomplished  by  hand  and  was  extremely  tedious.  If  these  techniques 
are  to  be  applied  in  the  future,  their  automation  is  clearly  desirable. 

The  use  of  one  cell  to  represent  a group  of  cells  with  varying 
characteristics  has  an  important  drawback;  as  the  variation  in  cell  char- 
acteristics increases,  the  percentage  of  cells  for  which  the  "representa- 
tive" cell  is  actually  representative  decreases.  Since  the  cells  in  the 
line  were  similar,  the  model  cells  compared  favorably  with  the  observed 
cells.  The  greater  variation  of  the  cells  in  the  cluster  resulted  in  a 
less  favorable  comparison.  The  indication  is  that  when  cells  have  substan- 
tially different  characteristics,  they  should  be  divided  into  sub-groups 
prior  to  idealization. 

The  results  of  this  study  suggest  that  the  simple  two-cylinder 
model,  when  appropriately  initialized,  can  simulate  many  of  the  character- 
istics of  observed  cells  remarkably  well. 

It  appears  that  although  horizontal  cell  dimensions  are  not  unim- 
portant, the  effect  of  their  variation  within  reasonable  ranges  is  limited 
to  minor  changes  in  cell  height  and  changes  of  several  decibels  in  reflec- 
tivity. The  variable  exerting  the  greatest  influence  on  the  development 
of  model  cells  was  the  relative  humidity  profile.  Furthermore,  the  results 
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pointed  up  the  inadequacy  of  using  relative  humidity  profiles  determined 
from  soundings  taken  in  the  open  atmosphere.  The  use  of  such  soundings 
is  based  on  the  unrealistic  assumption  that  the  cells  are  isolated  in  the 
larger-scale  environment.  Apparently,  the  effects  of  mesoscale  processes 
must  be  taken  into  account. 

One  way  to  include  the  effects  of  mesoscale  processes  is  to  take 
soundings  in  the  immediate  vicinity  of  the  cells,  preferably  in  the  weak 
precipitation  surrounding  the  updrafts.  This  would  be  the  best  way  to 
use  the  model  in  its  present  form,  but  it  is  highly  impractical. 

A second  method  would  be  to  utilize  assumptions  which  take  into 
account  the  effects  of  mesoscale  processes  on  the  humidity  profile  (this 
method  was  applied  in  this  study  when  the  open-atmosphere  relative  humi- 
dity soundings  were  altered).  The  assumptions  could  be  based  on  intuitive 
judgment,  crude  calculations  of  mesoscale  effects,  or  trial-and-error. 

In  any  case,  these  approaches  essentially  force  the  atmosphere  to  simulate 
the  model,  whereas  the  opposite  relationship  is  desired.  Furthermore,  it 
is  questionable  whether  anything  important  could  be  learned  about  the 
effects  of  the  cell  on  its  environment  if  that  environment  has  to  be  alter- 
ed to  produce  the  cell. 

The  best  approach  to  including  the  effects  of  mesoscale  processes 
in  the  model  would  be  to  include  the  cell  model  as  an  integral  part  of  a 
mesoscale  model.  This,  of  course,  is  a problem  of  much  larger  scope;  but, 
such  a model  would  provide,  in  addition  to  more  realistic  simulation  of 
cumulus  cells,  insights  into  cumulus-mesoscale  interactions. 

A further  refinement,  which  should  be  incorporated  into  the  larger- 
scale  model  discussed  in  the  previous  paragraph  and  could  be  incorporated 
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into  the  present  model,  would  be  the  addition  of  a realistic  modeling  of 
boundary-layer  effects.  There  is  reason  to  bel’’eve  that  much  of  the  air 
that  is  drawn  into  the  base  of  an  updraft  converges  from  the  boundary 
layer  surrounding  the  cell  (see,  e.g.  Zipser,  1978).  However,  in  the 
present  model  formulation,  the  air  that  enters  the  updraft  near  its  base 
is  entrained  from  the  outer  cylinder,  It  is  uniikely  that  this  air,  which 
has  subsided  in  the  outer  cylinder,  bears  much  resemblance  to  boundary- 
layer  air. 

The  addition  of  a boundary  layer  to  the  model  should  incorporate 
the  formulation  of  a less  artificial  initiating  impulse.  The  height  of 
the  impulse  as  well  as  its  temperature  and  humidity  should  be  geared  to 
boundary  layer  structure  and  conditions.  The  experiments  in  which  the 
final  relative  humidity  from  one  cell  was  used  to  initialize  a subsequent 
cell  illustrated  the  possible  contribution  to  the  humidification  of  the 
atmosphere  by  the  humidity  impulse.  In  experiments  of  this  nature,  in 
which  the  long-term  (relative  to  the  cumulus  time  scale)  effects  of  con- 
vection on  the  environment  are  investigated,  it  is  imperative  that  the 
most  realistic  humidity  perturbation  attainable  be  employed.  Further 
research  is  needed  to  determine  realistic  forms  and  magnitudes  of  the 
velocity  perturbation,  but  these  perturbations  should  be  based  on  observed 
vertical  velocity  fluctuations  in  the  boundary  layer,  reinforced,  perhaps, 
by  ascent  from  mesoscale  convergence. 
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